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Abstract 

In this paper we show how to relate European call and put op- 
tions on multiple assets to certain convex bodies called lift zonoids. 
Based on this, geometric properties can be translated into economic 
statements and vice versa. For instance, the European call-put parity 
corresponds to the central symmetry property, while the concept of 
dual markets can be explained by reflection with respect to a plane. It 
is known that the classical univariate log-normal model belongs to a 
large class of distributions with an extra property, analytically known 
as put-call symmetry. The geometric interpretation of this symmetry 
property motivates a natural multivariate extension. The financial 
meaning of this extension is explained, the asset price distributions 
that have this property are characterised and their further properties 
explored. It is also shown how to relate some multivariate asymmet- 
ric distributions to symmetric ones by a power transformation that 
is useful to adjust for carrying costs. A particular attention is de- 
voted to the case of asset prices driven by Levy processes. Based on 
this, semi-static hedging techniques for multiasset barrier options are 
suggested. 
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tion; semi-static hedging 

AMS Classifications: 60D05; 60E05; 60G51; 91B28; 91B70 



1 



1 Options and zonoids: an introduction 



The stop-loss transformation from mathematical insurance theory associates 
a random variable ( with its stop-loss E(£ — k)+ at level k. Here E denotes 
the expectation and x + = max(z, 0) for a real number x, where k is usually 
interpreted as the excess (part of claim which is not paid by the insurer). 

From the mathematical finance viewpoint, the stop-loss transformation 
is identified as the expected payoff from a European call option for ( = 
St = Soe rT r) being the price of a (say non-dividend paying) asset at the 
maturity time T, where So is the spot price and e rT r] is the factor by which 
the price changes, r is the (constant) risk-free interest rate and 77 is an almost 
surely positive random variable. In arbitrage-free and complete markets the 
expectation can be taken with respect to the unique equivalent martingale 
measure, so that the expected value of the discounted payoff becomes the 
call price. If the underlying probability measure is a martingale measure, 
then Ei] = 1 and the discounted price process S t e~ rt , t G [0, T], becomes a 
martingale. Unless indicated by a different subscript, all expectations in this 
paper are understood with respect to the probability measure Q, which is 
not necessarily a martingale measure. In this paper we do not address the 
choice of a martingale measure in incomplete markets. 

The expected call payoff E(SV — k) + can be considered a function of the 
bivariate vector (k,F), where 



is the theoretical forward price on the same asset with the same maturity, i.e. 
E(SV — k) + = E(Fi] — k) + . Deterministic dividends or income until maturity 
can be incorporated in the forward price, e.g. by setting F = Soe^ r ~ q ^ T in 
case of a continuous dividend yield q. 

When working with n assets, we write i] for an n-dimensional random 
vector (771, .. . , r] n ) such that the price Sn of the ith asset at time T equals 
Fifji with Fi being the corresponding forward price. We denote this shortly 
as 

S T = F o rj = (Firji, . . . , F n rj n ) . (1.1) 

In order to relate the expected payoffs to certain convex sets we need the 
following basic concept from convex geometry. 
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Definition 1.1 (see Sec. 1.7). The support function of a nonempty 
convex compact set K in the n-dimensional Euclidean space R™ is defined by 

h K {u) = sup{(x, u) : x G K} , net", 

where (x,u) is the scalar product in R n . 

For instance, if K — [ - 3/ j IS the line segment in R n with end-points ±x, 
then hx(u) = \ (x,u)\; if K = [0, x], then 1ik{u) = (x, u) + ; if X is the triangle 
in R 2 with vertices (0, 0), (a, 0), and (0, b), then hx{u) = max(uia, u^, 0) for 
all u = (iti, m 2 ) G R 2 . 

A function g : R n i— > R is called sublinear if it is positively homogeneous 
{g{cx) = cg(x) for all c > and x G R n ) and subadditive (g(x + y) < 
g(x) + g(y) for all x,y G R ra ). It is well known in convex geometry that 
support functions are characterised by their sublinearity property and that 
there is a one-to-one correspondence between support functions and convex 
bodies, i.e. nonempty compact convex subsets of R n , see e.g. [36, Th. 1.7.1]. 

With each integrable n-dimensional random vector 77 = (j]x, . . . ,r] n ) it 
is possible to associate a (n + l)-dimensional convex body which uniquely 
describes the distribution of 77. For this, consider (n+ l)-dimensional random 
vector (1,77) obtained by concatenating 1 and rj or, in other words, by lifting 
7] with an extra coordinate being one. In the financial setting this extra 
coordinate represents a riskless bond. Because of the lifting, we number 
the coordinates of (n + l)-dimensional vectors as 0,1, ... ,n and write these 
vectors as (u , u) for u 6l and u G R n or as (m , U\, . . . , u n ). 

Let X be the random set being the line segment in R n+1 with end-points 
at the origin and (1,7/), see [26] for detail presentation of random sets theory. 
The support function of X is given by 

h x (u , u) = max(w + U1V1 H ^ u nVn, 0) = ((u , u), (1, 77)) + 

for (uo, u) G R ra+1 . The integrability of 77 implies that hx{uo, u) is integrable. 
The expected support function E/ix is sublinear and so is the support func- 
tion of a convex body EX called the (Aumann) expectation of X, see [2"6"| 
Sec. 2.1]. For our choice of X, the set EX is called the lift zonoid of r\ and 
denoted by Z n . It is known that Z n determines uniquely the distribution of 
7], see [29| Th. 2.21]. Note that the zonoid of 77 appears from a similar (non- 
lifted) construction as the expectation of the random segment that joins the 
origin and 77, see [29| Th. 2.8]. 
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In the univariate setting we assume that n = 1 and rj = St/ F is a positive 
random variable. Then 

{u + uxEi] w , ui>0, 
u ,ui<0, (1.2) 

E(ii + Uii]) + otherwise 

for all (uq, U\) G M 2 . Figure [T] shows the lift zonoid of rj for various volatilities 
(0.25,0.5,0.75) in the log-normal case with E77 = 1 calculated for T = 1. 
The higher the volatility the larger (thicker) becomes the lift zonoid. The 
upper and lower boundaries of lift zonoids are the so-called generalised Lorenz 
curves, which can be easily parametrised, see [29], PP- 43 and 44]. 




Figure 1: The segment X = [(0, 0), (1, 77)] and the lift zonoid Z v for log- 
normal 77 with volatilities a = 0.25,0.5,0.75, drift /x = — \a 2 and maturity 
T = 1. 

Since the lift zonoid uniquely determines the distribution of random vec- 
tor rj = (rji, . . . ,r) n ), it also determines prices of all payoffs associated with 
St = (Sti, • • ■ , Stu), assuming that the underlying probability measure is a 
martingale measure. For instance, in the univariate case hz v {—k,F) (resp. 
hz v (k, —F)) is the non-discounted price of a European call (resp. put) option 
with strike k. A similar interpretation holds for basket options. The support 
function determines uniquely the lift zonoid, so that the prices of European 
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vanilla options (basket calls and puts) determine uniquely the distribution 
of the assets and so prices of all other European options. In the univariate 
case this fact was noticed by Ross [31], Breeden and Litzenberger [5], while 
Carr and Madan [12] presented an explicit decomposition of general smooth 
payoff functions as integrals of vanilla options and riskless bonds. In view of 
the positive homogeneity of support functions and central symmetry of lift 
zonoids, see [221 Prop. 2.15], it suffices in fact to have all call prices with 
parameter vectors (—k, u) with norm one in R n+1 , where k > and u stands 
for the vector containing the products of the weights and forward prices of 
the corresponding assets in the components. Alternatively, it suffices to have 
all call prices for any fixed k > and any u G W 1 . 

As just mentioned it is well known that the lift zonoid is centrally sym- 
metric. Section [2] begins by showing that the central symmetry property of 
lift zonoids is a geometric interpretation of the call-put parity for European 
options. 

The main question addressed in this paper concerns further symmetry 
properties of lift zonoids, their probabilistic characterisation and financial 
implications. Section [2] continues to show that in the univariate case (where 
lift zonoids are planar sets) the reflection at the line bisecting the first quad- 
rant corresponds to the dual market transition at maturity. In view of this, 
random variables that lead to line symmetric lift zonoids are called self-dual. 
This property has an immediate financial interpretation as Bates' rule [I] or 
the put-call symmetry [HI [TU] . For instance, the lift zonoids of the log-normal 
distribution in the risk-neutral setting (see Figure [1]) are line symmetric, 
which implies the put-call symmetry (or Bates' rule) for the Black-Scholes 
economy. Section [2] then shows how to translate the geometric symmetry 
property into symmetry relationships for general integrable payoffs and, in 
particular, for various binary and gap options. 

Section [3] highlights relationships between vanilla options and options on 
the maximum of the asset price and a strike. This leads to a concept of 
lift max-zonoids, which are particularly useful to describe options involving 
maxima of possibly weighted assets. This section also deals with a symmetry 
property of lift max-zonoids and shows how to relate option prices to certain 
norms on IR 2 yielding a relationship to the extreme values theory. 

Section H] characterises random vectors that possess symmetry properties 
generalising the classical put-call symmetry for basket options and options 
on the maximum of several assets. The symmetry (or self-duality) is under- 
stood with respect to each particular asset or for all assets simultaneously. 
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Relationships between the self-duality property and the swap-invariance in 
Margrabe type options have been studied in |28j. In currency markets, the 
self-duality results can be interpreted with respect to real existing markets, 
yielding the basis for further applications, see [35J . The overall symmetry im- 
plies that expected payoffs from basket options are symmetric with respect to 
the weights of particular assets and the strike price. We show that symmetries 
for some vanilla type options (like Bates' rule in the univariate case) imply a 
certain symmetry for every integrable payoff function. After discussing some 
fundamental results, we characterise the multivariate log-infinitely divisible 
distributions, exhibiting the multivariate put-call symmetry. The new effect 
in the multivariate setting is that independence of asset prices prevents them 
from being jointly self-dual. In other words, symmetry properties for several 
assets enforce certain dependency structure between them, which is explored 
in this paper. 

In order to extend the application range of the self-duality property and 
also in view of incorporating the carrying costs, we then define quasi-self- 
dual random vectors and characterise their distributions. These random 
vectors become self-dual if their components are normalised by constants 
(representing carrying costs) and raised to a certain power. The related 
power transformation was used in [TTI Sec. 6.2] in the one- dimensional case. 
Here we establish an explicit relationship between carrying costs and the 
required power of transformation for rather general price models based on 
Levy processes. 

These results are then used in Section [5] to obtain several new results 
for self-dual random variables thereby complementing the results from [TT] . 
In particular, self-dual random variables have been characterised in terms 
of their distribution functions; it is shown that self-dual random variables 
always have non-negative skewness and several examples of self-dual random 
variables are given. 

As in the univariate case also in the multiasset case there are various 
applications of symmetry results. First, symmetry results may be used for 
validating models or analysing market data, e.g. similarly as in [4] and [T7] 
in the univariate case. Furthermore, they could be used for deriving certain 
investment strategies, see e.g. Section |6~T1 The probably most important ap- 
plication will potentially be found in the area of hedging, especially in devel- 
oping semi-static replicating strategies of multiasset barrier and possibly also 
more complicated path-dependent contracts. Following Carr and Lee [TT] . 
semi-static hedging is the replication of contracts by trading European-style 
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claims at no more than two times after inception. As far as the relevance of 
this application is concerned we should mention that there has been a liquid 
market in structured products, particularly in Europe. At the moment the 
majority of the trades is still over-the-counter, but more and more trades are 
also organised at exchanges, especially at the quite new European exchange 
for structured products Scoach. Structured products often involve equity in- 
dices, sometimes several purpose-built shares, and quite often have barriers. 
Hence, developing robust hedging strategies for multi-asset path-dependent 
products seems to be of a certain importance. In the univariate case, Carr 
et al. [HI [9j [101 E] and several other authors (see e.g. [21 [lj [30]) developed a 
machinery for replicating barrier contracts having fundamental relevance for 
other path-dependent contracts. 

Section [6] contains first applications of the multivariate symmetry proper- 
ties, especially for hedging complex barrier options, thereby extending results 
from [HI [9j [10] and [TTJ for some multiasset options. The development of a 
more general multivariate semi-static hedging machinery is left for future 
research. 

2 Symmetries of lift zonoids and financial re- 
lations for a single asset case 

2.1 Parities 

We write c(k, F) for the price of the European call option with strike k on 
the asset with forward price F. Furthermore, let p(k, F) denote the price of 
the equally specified put. The maturity time T is supposed to be the same 
for all instruments. 

One of the most basic relationships between options in arbitrage-free mar- 
kets is the European call-put parity. In case of deterministic dividends, this 
parity can be expressed by 

c (k,F) =e~ rT {F -k)+p{k,F). (2.1) 

Recall that rj is defined by St = Frj, where St is the asset price at 
maturity and F is the forward price. The lift zonoid Z v of rj is centrally 
symmetric around | (1, E77), see [291 Prop. 2.15]. If the expectation is taken 
with respect to a martingale measure, then E?7 = 1, whence Z vo = Z v — (| , | ) 
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is origin symmetric, so that h Zrjo (u) = hz V}0 (~ u ) f° r an u e Interpreting 
the values of the support function of Z v as non-discounted call and put prices, 
this symmetry yields that 

e rT c(k,F) = h Zv (-k,F) = h Zri o+ (i ,!)(—&, F) 

= h z ^(k,-F)-h+ 1 -F 

= h Zv:0 (k,-F) + h-±F-k + F 

= h Zv (k, -F)-k + F = e rT p(k, F) + F-k, 

i.e. we arrive at the classical European call-put parity. By defining appropri- 
ate multidimensional lift zonoids and using their point symmetry, the above 
proof can easily be extended to the call-put parity for Asian options with 
arithmetic mean and to the European call-put parity for basket options. 



1.2-1 




") 0.2 04 0.6 0.8 1 1.2 



Figure 2: An approximation of the payoff set A for the Black-Scholes economy 
with volatility a = 0.5, interest rate r = 0.12, dividend yield q = and 
maturity T = 1. 

It is also easy to explain geometrically why the parities do not hold for 
American options. Let C(k,F) (resp. P(k ) F)) be the price of an American 
call (resp. put) with strike k on the asset with forward price F. Since the 
functions C and P are sublinear, it is possible to define convex body A with 



8 



support function h A (-k,F) = e rT C(k,F) and h A (k, -F) = e rT P(k,F), for 
k, F > 0. This convex body (that we call a payoff' set) determines the values 
of American vanilla options and is, up to some rare exceptions, not centrally 
symmetric, see Figure [2j 



2.2 Duality 

Recall that rj defined from St = Fr] is almost surely positive. If n is dis- 
tributed according to a martingale measure Q, then a new probability mea- 
sure Q can be defined from 

dQ _ 
dQ~ V ' 

Since rj is usually represented as c Ht for a semimartingale H t , < t < T, 
and — H t is called the dual to H t , the random variable rj = ii" 1 is said to be 
the dual of rj. The dual lift zonoid is defined as the Q-expectation of the 
segment that joins the origin and {l,fj). 

Lemma 2.1. If Z v is the lift zonoid generated by almost surely positive 
random variable n with En = 1, then 

7—7 

where Z v denotes the reflection of Z v with respect to the line {(u , Ui) 6 I 2 : 

Proof. For (u , «i) G M 2 

h Zfj (u ,Ui) = Eq(w + Uifj) + = E[(u + u^ 1 )^] 
= h Zv (ui,u ) = hz^u^ui) . 

noticing that the support function of Z n is obtained from the support function 
of Z n by swapping the coordinates. □ 

Since Z n is centrally symmetric with respect to (| , \ ), the set Z^ can also 
be obtained by reflecting with respect to the line {(uo, ui) : u\ = 1 — uo}. 

Lemma 12.11 relates the symmetry property of lift zonoids to the dual- 
ity principle in option pricing at maturity. This principle traces its roots 
to observations by Merton [22] > Grabbe [2D], McDonald and Schroder 
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Bates [3J, and Carr [B] and has been studied extensively over the recent 
years, e.g. Carr and Chesney [7J and Detemple [H] discuss American version 
of duality. For a detailed presentation of the duality principle in a general 
exponential semimartingale setting and for its various applications see Eber- 
lein et al. [15] and the literature cited therein. The multiasset case has been 
studied in Eberlein et al. [TB] . 

Consider a European call option valued at c(k, F, Q) with strike k and 
maturity T on a share represented in a risk-neutral world by a Q-price process 
St = <So e ? ^ t? 7t = Ftf]ti i-e. the share is traded in a market with deterministic 
risk-free interest rate r and attracts deterministic dividend-yield q, while r]t 
is a Q-martingale. In the most common setting this call option is related to a 
dual European put p(S , F, Q) with strike S and the same maturity T on a 
share represented by the dual Q-price process St = ke^ q ~ r ^ l f\t = F t f\ u where fj t 
is a Q-martingale. In other words, the dual put is written on another (dual) 
share traded in the dual market with risk-free interest rate q assuming that 
this dual share attracts dividend-yield r. By Lemma I2TT1 we get the following 
geometric proof and interpretation of the European call-put duality, 

p(S , F, Q) = e^ T EQ(5 - Ff,) + = e~ qT h Zfi (S , —F) 
= e-« T h 2v (S ,-F) 

= e -« T h Zv (-F, S ) = e-* T E Q (S oV - F) + 

= e- rT E Q (S , e^ T r ? - ke^ T+ ^ T ) + = c(k, F, Q) . 

By the classical call-put parity, a call-call duality can be derived in a straight- 
forward way. The duality for American options (see (TUTU [37]) can be inter- 
preted geometrically by reflecting the payoff set A (see Figure [2]) at the line 
bisecting the first quadrant. 

2.3 Symmetries for vanilla options 

If Z v is symmetric with respect to the line {(uq,ui) : Mo = ui} bisecting the 
first quadrant, i.e. 

Z-q = Zrj > (2-2) 

the duality relates out-of-t he-money to certain in-the-money calls in the same 
market. We call a positive random variable r\ self-dual if (12.21) holds. This 
symmetry property (12. 2\i clearly depends on the probability measure used to 
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define the expectation. It implies that E77 = 1, i.e. in the one period setting 
each probability measure that makes n self-dual is a martingale measure. 

Theorem 2.2. Assume that the asset price at maturity of an asset with 
deterministic dividend payments is St = Frj with self-dual n. Then 

e rT c{k, F) + k = e rT c{F, k) + F . (2.3) 

Proof. Noticing (I2.2p . we have 

e rT c(k, F) = h Zv (-k, F) = h^F, —k) = h Zv {F, -k) = e rT p(F, k) . 

The proof is finished by applying the put-call parity. □ 

Since c(k, F) is positively homogeneous in both arguments, i.e. c(tk, tF) = 
tc(k, F) for t > 0, the symmetry relation (I2.3P can also be written as 

c(m) + e~ rT 
c(m _1 ) + e~ rT 

where m = F/k determines the moneyness of the call and c(m) = c(l,m). 
Corollary 2.3. Assuming A2.2\) . the following European symmetries hold 

p(k,F) = c(F,k), (2.4) 
e rT p(k, F) + F = k + e rT p(F, k) . (2.5) 

Proof. By ([2l| and (fUD we obtain 

c(F, k) = c{k, F) - Fe~ rT + ke~ rT = p{k, F) . 

By combining the left-hand side of (12. 3p with (12.41) for the reversed order of 
k and F and the right-hand side of (12.31) with (12. 4p we arrive at (I2.5p . □ 

The above relations are known in the literature as put-call symmetry, see 
e.g. [1], [6], [10], and more recently [T71 [18] for log- infinitely-divisible mod- 
els, which are further discussed in Section 14.31 Further recent developments 
are presented in [11] . There are various applications of the put-call symme- 
try, especially in connection with hedging exotic options, see [TOT [TT] and 
Section [HI 
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2.4 General symmetry 

The following result obtained in [TTJ Th. 2.2] (without use of lift zonoids) 
generalises the self-duality to a wide range of payoff functions. 

Theorem 2.4. An integrable random variable n is self-dual if and only if for 
any payoff function f : M + i— > R such that E\f(Fr])\ < oo for F > 0, 

Ef(F V ) = E[/(Fr / - 1 )r ? ] . (2.6) 

Proof. Changing measure from Q to its dual Q, we arrive at 

E Q f(F V ) = Eq^F^tT 1 ] = E^fiFfj- 1 )^} = E^fiFn- 1 ^] . 

Since lift zonoids uniquely determine distributions of random variables, the 
last equality holds by Lemma [27T1 in view of the symmetry assumption (|2.2|) . 
Conversely, if (12.61) holds for any integrable payoff-function, then it holds for 
any vanilla options, i.e. hz v (—k,F) = h Zr] (F,—k) = (—k,F) for every 
k > 0, implying (Q. □ 

Denote by BC(k c , F) and GC(k c , F) the arbitrage-free values of a binary 
call and gap call with maturity T and strike k c , i.e. the European derivatives 
with payoffs ls T> k c and St^St^- Furthermore, BP(k p ,F) and GP(k p ,F) 
denote the arbitrage-free values of the equally specified puts, i.e. the Euro- 
pean derivatives with payoffs l ST<kp and SV Is T< fc p . Theorem 12.41 yields the 
following result, which is equivalent to [11] Cor. 2.9] being a generalisation 
of the binary put-call symmetry from [10] . 

Corollary 2.5. Under the assumptions of Theorem \2.4\ the following rela- 
tionships hold: 

V / ^BC(A; C , F) = -^= GP(fc p , F), ^k~ p BP(k p , F) = GC(fc c , F) , 

where the forward price F equals the geometric mean of the binary (resp. 
gap) call strike k c and the gap (resp. binary) put strike k p , i.e. F = ^k c k p . 

Proof. If ^Jk~k p = F, then (JES) yields that 

e rT BC(k c , F) = E[l F71>kc ] = E[ V l Fv -i >kc ] = E f ' 



n — : 1 / — 

ll- I- \/k c k p >k c r) 

■ y 

1 rT 

-.E[S T l kp>ST ] = -4= GP(A; p , F) . 



\fk c kp \/k c k p 
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The proof of the second identity is similar. 



□ 



Thus, symmetry properties for particular options (vanilla, binary/gap, 
straddles, etc.) are nothing but writing down Equation (12. 6p for special 
payoff functions. 




Figure 3: Symmetries of Z n for self-dual r\ and their financial interpretations. 



The random variable r] has no atoms if and only if the support function 
hz is continuously different iable on R 2 \ {0}, so that Z v is strictly convex, 
see [36, Sec. 2.5]. Then the unique point on the boundary of Z n (without 
the points (0,0), (1, 1)) at which (uo,Ui) G R 2 \ {0} is the outward normal 
vector to Z r] is given by 

gradHK« 1 )=(^^ L ) • 

Assuming that the underlying probability measure Q is a martingale measure 
and calculating the gradient of the support function hz at (m , «i) = (—A;, F) 
with k, F > yield another parametrisation of the upper boundary of Z v as 

grad h Zv (-k, F) = e rT (BC(fc, F), GC ^ F A j k > . 
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Analogously, the lower boundary is parametrised by 




k > 0. 



Hence, in non-atomic cases the boundaries of the lift zonoid Z v can be 
parametrised by non-discounted arbitrage- free values of binary and normalised 
gap options. Thus, the distribution of r] is also reflected in these pairs of op- 
tions. 

Figure [3] interprets financial parities and symmetries geometrically for the 
lift-zonoid Z v . By comparing the coordinates of the points A, B, C, D we get 
relations between binary and gap options. By comparing the related values 
of the support function hz„ we arrive at the put-call parity and symmetry 
for vanilla options. Combining A with C yields parities between certain in- 
the-money puts (vanilla, binary, and gap) and the related out-of-the-money 
calls. Connecting the points B and D yields the same parity results, except 
that B represents certain in-the-money calls and D the related out-of-the- 
money puts. Comparing C with D yields out-of-the-money put call symmetry 
results, while linking A with B leads to the same results for in-the-money 
options. Finally combining B with C (resp. A with D) results in the call-call 
(resp. put-put) symmetry. 

3 Options on maximum and lift max-zonoids 

Since the call payoff can be written as 



the expected call payoff can be related to another convex compact subset M v 
of M 2 that has the support function 



The set M v is defined as the Aumann expectation of the random triangle with 
vertices at the origin, (1,0), and (0, if). Because the financial quantities are 
non-negative, we often restrict the support function onto the first quadrant 
W 2 + . Then it is possible to write (13.21) as 



(St — k) + = max(S , T, k) — k , St, k > , 



(3.1) 



h Mv (k, F) = E max(fc, Frj, 0) , (k, F) G R 2 . 



(3.2) 



h Mv (k,F) =Em&x(k,Fri) , (k, F) G M 2 



(3.3) 
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If i] = . . . , n n ) is a random vector in W] = [0, oo) ra , a similar to (13.21) 
construction leads to the set M, n called the lift max-zonoid of rj and defined 
as the Aumann expectation of the random crosspolytope with vertices at the 
origin and the unit basis vectors eo, e%, . . . , e n in R n+1 scaled respectively by 
Mi, '■■iVn, i-e. 

hM n (uo,ui, ...,u n ) = Emax(0,M ,ui?7i, . . .,u n r) n ) , (u ,ui, ...,«„) G R n+1 . 

Max-zonoids have been introduced in [27] in view of their use in extreme 
values theory. If Er/ = (1, . . . , 1), then M v is a convex compact subset of 
the unit cube [0, that contains the origin and all unit basis vectors. If 
n = 1, then each such set is a lift max-zonoid of some random variable, while 
this no longer holds for two and more assets, see [271 Th. 2}. 

Theorem 3.1. The lift max-zonoid M v of an integrable random vector rj G 
determines uniquely the distribution of r). 

Proof. The support function hM v (uo, Ui, . . . , u n ) for w , u%, . . . ,u n > can 
be written as Emax(« ,() for ( = max(wi?7i, . . . ,u n n n ) and so determines 
uniquely the distribution of (. The cumulative distribution function of ( is 
given by 



and so determines uniquely the joint cumulative distribution function of 



If the underlying probability measure is a martingale measure, Theo- 
rem 13.11 implies that prices of options on the maxima of weighted assets 
determine the joint distribution of the risky assets and so prices of all other 
payoffs. In view of the positive homogeneity of support functions it suffices 
that the expected values are known for parameter vectors (uq, u) with norm 
one in M n+1 , m > and u with strictly positive coordinates. Alternatively, 
it suffices to know the expected values for fixed Uq > and u with strictly 
positive coordinates. 

The remainder of this section deals with the single asset case. Equation 
(13.11) suggests that in this case, the lift max-zonoid is closely related to the 
lift zonoid Z n of a random variable rj. 



F c {t) = P{ Vl < 




□ 
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Lemma 3.2. If M v is the lift max-zonoid generated by a non-negative inte- 
grable random variable rj, then 

M, = conv{(0,0)U(Z; + (l,0))}, 

where Z' v is the reflection of Z v with respect to the line {(uq, u\) : u = 0} 
and conv denotes the convex hull, see Figured 

Proof. For (u ,Ui) G IR 2 the support function of the set in the right-hand 
side is given by 

max(0, h z ^(u , «i) + u ) = max(0, h Zr) (-u , u x ) + u ) 

= max(0, E(mi?7 — u o)+ + u o) ■ 

By checking all possible signs of uq and u\ it is easy to see that this support 
function equals Emax(0, uq, uirf) = hM v (uo,ui). □ 




Figure 4: Relation between M v and Z v as well as the body M v for log-normal 
r) with mean one and volatility o = 0.5 calculated for T = 1. 

Lemma 13.21 implies that the duality transform from Section 12.21 amounts 
to the symmetry of M n with respect to the line bisecting the first quadrant. 
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Furthermore, r\ is self-dual if and only if M v is symmetric with respect to 
this line, i.e. 

E max(Fr], k) = E max(F, krj) , k, F > . (3.4) 

Indeed, since 

hz v (- u o,Ui) = E maxOiT?, w ) - w , 
frz„(ui, -«o) = Emax(«i,M ?/) - tt , 

for Wo,«i > 0, the symmetry property (12.21) is equivalent to (13.41) . 

The Euclidean space IR 2 can be equipped with various norms. Each norm 
on M 2 can be described as the support function of a centrally symmetric (with 
respect to the origin) convex body that contains the origin in its interior. 
Although M n in Lemma 13.21 is a subset of IR+ and so does not contain the 
origin in its interior, it is possible to use hu to define the norm on the whole 
plane as 

\Mv = h M n {\x\) = h^{x) , 

where |x| is the vector composed of the absolute values of the components 
of x G M 2 and M v is obtained as the union of symmetrical transforms of 
M v with respect to the coordinate lines, see Figure HI For instance in the 
martingale setting the call price satisfies 

e rT c(k,F) = \\x\\ v - k, x = (k,F)eM. 2 + . 

Conversely, each norm on determines uniquely the distribution of an 
integrable non-negative random variable. 

Note that rj is self-dual if and only if the norm || • ||^ is symmetric, i.e. 
||(tto,«i)||,, = ||(iti,«o)||»7 for a11 ( u o,Mi) e R%. 

Example 3.3. Consider the £ p -norm on M 2 , which is clearly symmetric. Eval- 
uating the £p-norm of (t, 1), we arrive at 

xp r ,(x)dx, t>0, 

assuming that rj is absolutely continuous with density p v . Differentiating 
with respect to t yields that 

P( v <t) =t p - 1 {t p + l) 1/p - 1 . 
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Thus, 77 has the density 



l/p-2 



i > 0, 



which is shown to imply the self-duality of 77, see Corollary 15.2( a). 

We conclude this section by stating a relation between norms and extreme 
values. It is known [191 EZ] that each norm || • || on corresponds to a 
bivariate max-stable random vector (£1,^2) with unit Frechet marginals, i.e. 



An important norm on M 2 related to the Black-Scholes formula and the theory 
of extreme values is mentioned in Example 15.31 

4 Multivariate symmetry 

4.1 Characterisation of distributions with symmetry 
properties 

It has been shown in Section [2] that geometry of the lift zonoid for a single 
asset price has a fundamental financial importance. The lifting operation 
amounts to adding an extra coordinate to the asset price or prices and so 
increases the dimension by one. For instance, the lift zonoid associated with 
the integrable price of a single asset is a subset of the plane, which is always 
centrally symmetric convex and compact. It is well known [36] that all cen- 
trally symmetric planar compact convex sets are zonoids, while this is not 
the case in dimension 3 and more. This fact already suggests an important 
dimensional effect that appears when dealing with more than one asset. 

For the multiasset case the direct relationship between lift zonoids and 
lift max-zonoids is also lost. Indeed, the maximum of two numbers can be 
related to the stop- loss transform as max(a, b) = (a — b) + + b, while this no 
longer holds for the maximum of three numbers. The family of multivariate 
symmetries is also considerably richer than in the planar case. 

Let St = (Sti, • • • , Sm) and integrable r/ = (7/1, ... , r] n ) be as defined 
in (11. ip . Assume that all coordinates of r\ are positive, i.e. 77 6 E n = (0, oo) n 
and rj = e^ for a random vector £ = (£1, ...,£«), where the exponential 
function is applied coordinatewisely. For simplicity of notation, we do not 
write time T as a subscript of rj and incorporate the forward prices Fj, 



P(6 < u 



-1 
1 j 



£2 < u 2 l ) = exp{-||(wi,w 2 )||}, (ux,u 2 ) e 
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j = 1, . . . ,n, into payoff functions, i.e. payoffs will be real-valued functions 
of Tj. 

In the sequel two specific payoff functions are of particular importance, 
namely 

n 

f b (u , ui,..., u n ) = ( ujrji + up) , uq, Ui, . . . ,u n G K , 

1=1 + 

for a European basket option and 

n 

f m {u , Ui, ...,«„)= liQ V Y tljT/i , Uo, «i, . . . , u n > , 

Z=l 

for a European derivative on the maximum of n weighted risky assets together 
with a riskless bond, where V denotes the maximum operation. Despite 
the fact that the payoff functions /& and f m depend on 77, we stress their 
dependence on the coefficients, since it is crucial for symmetry properties. 

By ( 13. ip . call and put options on the maximum of several assets can be 
written by means of the payoff function f m , e.g. 

n 

(\Juirji-kj = f m (k,ui, . . . ,u n ) - k , k>0. 
1=1 + 

In view of Theorem 13.11 prices of these options uniquely characterise the 
distribution of an integrable random vector rj e . 

If X is the segment that joins the origin in M n+1 and (1,77), then /& 
becomes the support function of X, so that the expected payoff is the support 
function hz v (u , ui, . . . , u n ) of the lift zonoid Z v , i.e. 

E/&(uo,«i, •••a) = h Zrj {u Q ,u u ...,u n ). 

Similarly, for it , iti, . . . , u n > 0, the expectation of f m becomes the support 
function of the lift max-zonoid M v . 

Fix an arbitrary asset number % G {l,...,n} and assume that Q is a 
probability measure that makes rj integrable. Recall that E without subscript 
denotes the expectation with respect to Q, otherwise the subscript is used to 
indicate the relevant probability measure. Since 7/ 1 / 2 = [i]^ 2 , • • • , f]n) — eh^ 
is integrable, we can define new probability measures Q* and 8 l by 

dQ l rji dS l ea& 

~d~Q ~ E^' dQ ~ Ee^ 1 ' 
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Hence, £ l is the Esscher (exponential) transform of Q with parameter ~ ej, 
where is the ith standard basis vector in R n , see [32] and [Ml Ex. 7.3] for 
the Esscher transform in the context of multivariate Levy processes. Since 

we see that Q is the Esscher transform of £ l with 

parameter — 

For simplicity of notation, define families of functions x, : E n i— > E n and 
linear mappings K; t : W 1 > IR n acting as 

X\ X{—\ 1 Xi-\-\ x n 

Xj I i • • • i i i i • • • i 



K^x (^i x^, • • • , Xi—\ x^, x^, X{-^-\ x^, . . . , x n x^) (4.1) 

for i = l,...,n. The linear mapping Ki can be represented by the ma- 
trix Ki = (/c^ m )^ =1 with ku = 1 for all I ^ i, k u = — 1 for / = 1, . . . ,n 
with all remaining entries being 0. Note that Xj and Ki are self-inverse, i.e. 
Xj(xj(x)) = x and KiKiX = x, and that the ith coordinate of KiX is — x«. 
The transpose of Ki is denoted by Kj . In the following we consider vectors 
as rows or columns depending on the situation. 

The permutation of the zero-coordinate with the ith coordinate of a vector 
(u , u) G M. n+1 is denoted by 

TTi(u , U) = (U i} Ml, ... , Ui-i, M , ■■■,%) for i = 1, . . . ,71. 

If L? C M n+1 , then 7Ti(B) is the reflection of B at the hyperplane {(mo,m) G 
E n+1 : W j = m }- 

Finally, (resp. <y2| l ) denotes the characteristic function of the random 
vector £ under the probability measure Q (resp. £ l ). 

Univariate versions of the statements (i), (iii), (vi), and (vii) of the fol- 
lowing theorem are already known from [TTj Th. 2.2, Cor. 2.5]. 

Theorem 4.1. Let i] = be an n-dimensional Q-integrable random vector 
with positive components and let i be a fixed number from {1, . . . , n}. The 
following conditions are equivalent. 

(i) For all u ER and u G W 1 , 

Ef b (u ,u) = Ef b (iri(u ,u)) . 

(ii) For all u > and u G , 

Ef m (u Q ,u) = Ef m (-Ki(u ,u)) . 
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(Hi) For any payoff function f : E n i— > R such that E|/(?y)| < oo we have 



Ef{ v ) = E[/(*fo))ifc] . 

(iv) The lift zonoid Z v oft] satisfies iri(Z v ) = Z v . 

(v) The lift max-zonoid M v of rj satisfies ni(M v ) = M v . 

(vi) The distribution of 77 under Q is identical to the distribution of ff = 
Kiiyj) under Q\ 

(vii) The distributions of £ and under E % coincide. 

(viii) For every u G W 1 , 

yf{u) = yf{Kju) 

or, equivalently, 

ff (u - ~ ie^J = <ff (Kju - i te^j , 

where t = is the imaginary unit and e{ is the ith standard basis 
vector in IR n . 

Since (vi) corresponds to the duality transform in the univariate setting 
(see Section 12. 2p . we say that 77 satisfying one of the above conditions is 
self- dual with respect to the ith numeraire and write shortly 77 G SDj. If 
77 is self-dual with respect to all numeraires i = 1, . . . ,n, we call rj jointly 
self- dual. 

Remark 4.2 (Relaxing of (i) and (ii)). In view of the positive homogeneity 
of payoff functions /& and f m it suffices to impose (i) and (ii) for parameter 
vectors (v,q,u) with norm one in M n+1 , or, in (i), for any fixed Uq 7^ and 
any u G M n . Condition (ii) can be assumed only for any fixed m > and u 
with strictly positive coordinates, i.e. for u G E n . 

Remark 4.3 (Martingale property). If rj G SDj, then (iii) applied to / identi- 
cally equal one (or symmetry conditions (iv), (v), (vi) for lift (max-) zonoids) 
imply that E?7j = 1, i.e. Q is a martingale measure for the ith component of 
77 in the one-period setting. However, Q does not need to be a martingale 
measure for other components of 7], quite differently from the univariate case 
[TT] . The martingale property for all components is ensured by requiring that 
77 is jointly self-dual. Otherwise it has to be imposed additionally, if needed. 
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Remark 4.4. If the forward prices are not included in the payoff function, 
then condition (iii) for the self-duality of v with respect to the ith numeraire 
can be equivalently expressed as 



E/(S 



TI, 



E 



S Tn ) = Ef(Forj) = E[f(F o^(rj)) n , 
StiF{ 



Fi) 2 S T (i+i)Fi 



S- 



Ti 



Ti 



Ti 



Ti 



Sti Fi 



by applying (iii) to f(rj) =f(Forj). 

Remark 4.5 (Conditioning in Theorem 14.11) . All conditions of Theorem 14.11 
can be written conditionally on a fixed event or conditionally on a cr-algebra. 
This has the following application for stochastic processes. Consider a fam- 
ily {v{t), t > 0} of random vectors being self-dual with respect to the zth 
numeraire. If r is a non-negative random variable, which is independent of 
{v(t), t > 0} then r](r) satisfies all statements in Theorem 14. II with the ex- 
pectations taken conditionally on the a-algebra generated by r, i.e. v( T ) is 
conditionally self-dual with respect to the ith numeraire. 

To prove Theorem 14.11 we need the following multivariate extension of 
the duality principle at maturity by reflection, see Lemma 12. 1 1 The dual lift 
zonoid (resp. lift max-zonoid MfA with respect to the zth numeraire Vi is 
defined for the random vector rf = (fj\, . . . ,fj l n ) = Xi{v) and the expectation 
with respect to Q\ 

Lemma 4.6. Let rj = (rjx, ... , t] n ) be an integrable random vector with Er^ = 
1 for a fixed i G {1, . . . , n}. Then 

Zfji = TTi(Z v ) and = •n i {M r) ) . 

Proof. For (uq,u) G M n+1 we have 



n 

h Zfji {u , u) = E Ql ( U W\ + 



u 



E 



i=i 

n 



E Q d > ui 1 \-u 



Em Ui 
ui— H h u ) Vi 



Vi Vi 



1=1, tyi 
n 

E ( uivi + Ui + u v?j 

1=1, l^i + 

hz v (n(uo,u)) = K^ Zri )(uo,u) . 
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The proof for lift max-zonoids is similar. □ 

Proof of Th. \4-l\ We will establish all equivalences in several steps. 

(i)=r»(iv)=^(vi) The definition of the lift zonoid, (i) (also implying Erji = 1), 
and Lemma [4.61 imply that for all (uq,u) G R n+1 

h Zv (u ,u) = Ef b (u ,u) = Ef b (TTi(u ,u)) 

= h Zv (TTi(u ,u)) = h ni{Zv) (u ,u) = h Zfji (u Q ,u) , 

so that Z n and coincide (see (iv)) as having identical support functions. 
Since the lift zonoid uniquely determines the distribution of an integrable 
random vector, this implies (vi). 

(vi)=r>(iii)=^(i) The definition of Q\ the self-inverse property of x« and (vi) 
yield that 

E Q /(i7) = Evimvr 1 } = E Ql [/(x,(^))^] = E Q [f(^( V )) Vl ] , 

so that (iii) holds. By applying (iii) for the payoff function /& of a basket 
option we arrive at (i). 

(iii)=^(ii)=^(v)=^(vi)=^(iii) By applying (iii) to the payoff function f m we 
obtain (ii). Now the definition of the lift max- zonoid, (ii), and Lemma [4.61 
yield that 

hM v ( u o,u) = h Mv (^i(u ,u)) = h ni (M v ){u ,u) = h Mf)i (u ,u) 

for every (u ,u) G E ra+1 , and thus, for every (u ,u) G M n+1 , i.e. (v) and (vi) 
hold, since the lift max- zonoid uniquely identifies the distribution, by Theo- 
rem [3J3 It is already shown that (vi) implies (iii). 

(iii)<^>(vii)<^>(viii) Since m; = E(es^) is finite, 

Eftn) = E/(e*) = m^/^e-H] , 
E[/(x i (r 7 ))7 7i ] = E[/(e**)e*] = mi E £i [f{e K *)e^] . 

Thus, (iii) yields that 

E £i [/(e*)e-^] = E fi [/(e^)e^] = E £i [f(e K *)e-^ K ^} (4.2) 

for any Q-integrable payoff function /. Recall that the ith coordinate of 
is — £j. Choosing f(x) = g{x)e^ x \ we see that the ^-expectations of g(e^) 
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and g(e Ki ^) coincide for all continuous functions g with bounded support, 
whence £ coincides in distribution with K£ under £ l . Conversely, if £ and 
Kit; share the same distribution under S % , then (14. 2[) holds and implies (iii). 
Furthermore, (vii) is equivalent to 

<pf (w) = ¥>£rf(u) = <pf {Kju) 

for every u G M n . Writing the characteristic functions as ^-expectations and 
referring to the change of measure, the latter condition is equivalent to 



E 


it. 
„ it 


= E 




Ee^ 1 - 





<pf ( u - \ * e i) = <pf (kJu - ^ le i 



so that 

1 \ 0/-.T 1 



for all ti6l" □ 

In view of Theorem I4.1( iv.v). examples of random vectors rj G SDj can 
be derived by constructing lift (max-) zonoids which are symmetric with 
respect to the hyperplane {(uo, • • • , u n ) G IR n+1 : uq = see Examples 13.31 
and 14.131 The following result can be helpful for such constructions. Its 
univariate version is stated in [371 Ex. 8]. 

Theorem 4.7. Consider an integrable random vector n G E n with distribu- 
tion Q. 

(a) If n is absolutely continuous with probability density p v , then r\ G SDj 
if and only if 

p v (x) = x^ n+2 ^ p v (xi(x)) for almost all x G E n , (4.3) 
equivalently, the density p^ of £ = log rj satisfies 

Pt{x) = e~ Xi p^(Kix) for almost all x G M™ . (4.4) 

(b) lfr\ is discrete, then r\ G SDj if and only ifQ(r) = Xi{x)) = XiQ(r] = x) 
for each atom x of 'rj. 
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Proof, (a) Condition (iii) of Theorem 14.11 can be written in the integral form 



as 



f(x)p v (x)dx = I f(x t (y))y i p v (y)dy= / f(x)-^ pi p r) (x i (x)) dx 



where the last equality is obtained by changing variables x = Xi{y) and 
noticing that = x. 

Consider the function f(x) = lxe[ai,6i]x-x[on,6 n ] f° r an Y parameters 
ax, . . . , a n , bi, . . . , b n G R. Differentiating both sides with respect to bx, . . . , b n 
and by using the dominated convergence, we get ( 14. 3ft almost everywhere. For 
the converse, write the right-hand side of (iii) as integral, refer to ( 14.31) and 
change variables. The equivalence between (14.31) and (14.41) can be seen by 
the classical density transformation, 
(b) In the discrete case (iii) can be written as 

f(x l )Q(v = x l ) = /(x^))x'Q(77 = x*) , (4.5) 
i i 

where the sum stretches over all atoms x l with x\ being the zth component 
of x l . Since ( 14. 5 j) also holds for f(x l ) = ll x=x i for any fixed atom x, we obtain 
that x* = Hi{x) is also an atom with 

Q(t7 = x )= x*Q(r) = x*) = — Q(?7 = x;(x)) . 

Xi 

For the converse we have for any (v,q, u) G IR n+1 

n 

h Zv (U , M ) = E ( M ° + E UmXl m) QiV = xl ) 



m=l 



n ^ 
^2{ U + ^2u m X l m ^ -jQ(T] = Hi(x 1 )) 

I m=l * 

n 

Y [ u i + u o x i + Yl UmX ™) QC 7 ? = ^*) 



I* m=l, m^i 

where x l * = Xj(x'), so that we obtain (i) of Theorem 14.11 □ 
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Now we give a result about the marginal distribution of 77, for the random 
vector n being self-dual with respect to this numeraire. 

Lemma 4.8. IfrjE SDj, then rji is a self-dual random variable. 

Proof. Choose in Theorem 14.1( h) vector u with all coordinates being zero 
apart from u and u, t . Then (ii) reads E(u V u^i) = E(uj V u r]i) for every 
u ,Ui > 0. This is exactly the symmetry condition (13.41) . □ 

4.2 Jointly self-dual random vectors 

Recall that random vector rj is called jointly self-dual if it is self-dual with re- 
spect to all numeraires. Since permutations of coordinate and an arbitrary 
i G {1, . . . ,n} generate by successive applications the transpositions of any 
two i, j G {0, 1, . . . , n}, the expected payoff functions /& and f m for jointly 
self-dual rj are invariant with respect to any permutation of their arguments, 



for each permutation i i— > 1^. In view of this, Theorem 14. II implies the follow- 
ing result. 

Theorem 4.9. Random vector 77 is jointly self-dual if and only if its lift 
(respectively lift max- ) zonoid Z v (respectively M v ) is symmetric with respect 
to each hyperplane {(uq, u±, . . . , u n ) G IR n+1 : u\ = Uj} for all i, j — 0, . . . , n, 

Corollary 4.10. If rj is jointly self-dual, then all its components are iden- 
tically distributed self-dual random variables with expectation one and n is 
exchangeable, i.e. its distribution does not change after any permutation of 
its coordinates. 

Proof. The components of 77 are self-dual by Lemma 14.81 and so have expec- 
tation one. If r] is jointly self-dual, then Theorem 14.91 yields that for every 
(u Q , u) G and any i, j = 1, . . . , n, 



e.g. 



E/bOo, Mi, ... , u n ) = Ef b (u io ,ui 



1 j 



■■■,Uln) 



(4.6) 



n 



n 




1=1, l^i 



1=1, 
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By setting u\ = for all / 7^ 0, % we arrive at 

E(u V UiT]i) = E(m V u^j) for every (u , Mi) G R+ . 

Thus, for any i,j = 1, . . . ,n the random variables rji and 77^ have the same 
lift max-zonoid. By Theorem 13.11 all coordinates of 77 share the same distri- 
bution. 

Theorem 14.91 yields that 77 = (771, ... , n n ) and . . . , r]i n ) obtained by 
any permutation of its coordinates share the same lift max-zonoid and thus 
the same distribution, i.e. 77 is exchangeable. □ 

It should be noted that the converse statement to Corollary 14.101 does 
not hold, i.e. the exchangeability of 77 does not imply joint self-duality. This 
is easily seen as a consequence of the following result, which says that any 
non-trivial random vector 77 with independent coordinates cannot be jointly 
self- dual. 

Theorem 4.11. Assume that n > 2. 

(a) If 77 G SDj and rji and rjj are independent for some j 7^ i, then rji equals 
1 almost surely. 

(b) If r] is a jointly self-dual random vector with independent coordinates, 
then all coordinates of 77 are deterministic and equal 1 almost surely. 

Proof. It suffices to prove only (a). By Theorem I4.1( ii) letting Ui = for 

E(u V Uifji V Ujrjj) = E(v,i V u rn V utfj) 
for all Uo,Ui,Uj G K.+. In particular, if = 0, then 

E(u V Uji]j) = E(u 77j V UjT]j) for all (u , Uj) G IR+ . 

Since 77^ is self-dual by Lemma l4~8| the conditioning on rjj yields that E(-u ^j V 
Ujfjj) = E(w V UjiiiTij). Hence, 

E(u V UjTjj) = E(uo V Uji]i7]j) for all (uo, Uj) G , 

whence r/j coincides in distribution with 77^77^ , see Theorem 13.11 If = log 77/ 
for I = i, j, then £ 3 - and £j + £j share the same distribution. Therefore, the 
characteristic function of £j identically equals one for some neighbourhood of 
the origin, whence = almost surely. □ 
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Remark 4.12 (Random vectors sampled from Levy processes). Assume that 
Ci, • • • , Cn are independent integrable random variables. Consider the vector 

V = (Cl> C1C2, C1C2C3, • • • , Cl • • • Cn) ■ 

This construction is important, since then £ = log?7 is a vector whose com- 
ponents form a random walk. However, 77 cannot be jointly self-dual as a 
vector, unless in the trivial deterministic case. Indeed, setting u ,Ui,u 2 > 
and u 3 = ■ ■ ■ = u n = 0, writing the expected payoff f m conditionally on ( 2 
and using the self-duality of (which follows from r\ G SDi) we see that 

E(w V u\f]x V U21I2) = E(«oCi V Mi V M2C2) 

is symmetric in uq,Ui,u 2 if r] is jointly self-dual. Thus, (Ci, C2) is a jointly 
self-dual vector with independent components, which is necessarily trivial by 
Theorem 14.111 An extension of this argument shows that Qy = ■ ■ ■ = ( n = 1 
almost surely. Therefore, it is not possible to obtain jointly self-dual random 
vectors by taking exponentials of the values of a Levy process at different 
time points. 

Example 4.13. The most obvious convex body in M n+1 being symmetric with 
respect to the hyperplanes u = u iy i = 1, . . . , n, is the closed unit ball 
Bi(0) of radius one centred at the origin. The value of the corresponding 
derivative equals the discounted magnitude of the weight vector. It is shown 
in [27] that M v = Bi(0) n is a max-zonoid. By Theorem M v is 
the lift max-zonoid of a jointly self-dual random vector rj, such that for all 
(u ,u)eR n + +1 

It p QQ TL 

hM v (uo,u) = \\(u ,u)\\ = B\u V \/uirnJ = u + / P (\J uitji > tj dt 

1=1 Ju ° 1=1 

+ rii-F V (±...±))dt, 

Jun V U l U n ' 



U 



where is the joint cumulative distribution function of rj. Using the ex- 
pression for the Euclidean norm ||(wo, w)||, differentiating with respect to all 
components and setting uq = 1 yield the following expression for the density 
of 77 

2 n T(n + - ) 

p v (u) = 1- , u = (tlx, . . . ,u n ) e E n , 

v^(i+Er=i«r 2 )^ +n nr=i«f 



28 



where T(-) denotes the Gamma function. It is easy to check that p v satis- 
fies ( 14. 3 p for every i = 1, . . . ,n. 

Example 4.14. Let Co>Cij---;Cn be i.i.d. self-dual random variables (their 
examples are provided in Section [5]). Define rji = (oQ for i = l,...,n. 
Conditioning on £i, • • • , Cn yields that 

EOo V UiTfr V • • • V U n 7] n ) = E[E(m V Co(wiCl V • • • V U„Cn)|Cl> Cn)] 

= E(m Co V UiCl V • • • V M n Cn) 

is symmetric in it , Mi, . . . , u n , i.e. r\ = (rji, ... , r^ n ) is jointly self-dual. Note 
that 7/1, . . . , Tj n are all self-dual random variables, but are no longer indepen- 
dent. In particular if the £'s are log- normally distributed with \i = — | and 
a — 1, then log 7/ is normally distributed with mean (— | , . . . , — | ) and the 
covariance matrix having diagonal elements one and all other | . We will 
return to this situation in Example 14.191 

4.3 Exponentially self-dual infinitely divisible random 
vectors 

A random vector £ has an infinitely divisible distribution if and only if £ = Li 
for a Levy process L t , t > 0, see [33]. In view of the widespread use of Levy 
models for derivative pricing we aim to characterise infinitely divisible ran- 
dom vectors £ = log r\ for r\ being self-dual with respect to the zth numeraire 
or all numeraires. If r] G SDj, then £ is said to be exponentially self-dual with 
respect to the ith numeraire and we write shortly £ G ESD.;. 

The Euclidean norm || • || is not invariant with respect to the transforma- 
tion x h-> KiX defined by ( 14. lft . For simplifying the formulation of the results 
we introduce the following norm on M™ 

Nll 2 = ^(lkll 2 + ll^ll 2 ), (4.7) 

where the number i G {1, . . . , n} is fixed in the sequel. It is easy to see that 
I'll is indeed a norm, which is equivalent to the Euclidean norm on MJ 1 . Since 
Ki is self- inverse, |||u||| = |||-Kjw||| for every u G M n . 

We use the following formulation of the Levy-Khintchine formula, see [33| 
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Ch. 2] , for the characteristic function of £ 

<pf(u) = Ee lM = exp ji( 7 , u) - ^ (u, Am) 

+ / (e^-l-t^x)!,^!)^)), (4.8) 

for u G M n , where A is a symmetric non-negative definite n x n matrix, 
7 G M. n is a constant vector and v is a measure on M. n (called the Levy 
measure) satisfying ^({0}) = and 

/ mm(\\x\\ 2 ,l)du(x) < oo . (4.9) 

Note that the latter condition can be equivalently written in the new norm 

Theorem 4.15. Let rj he an integrable random vector under probability 
measure Q such that £ = log rj is infinitely divisible under Q. Then £ G ESDj 
if an only if for the generating triplet (A, u, 7) the following three conditions 
hold. 

(1 ) The matrix A = (%)^ =1 satisfies = = \ an for all j — 1, ... , n, 
j ^ i- 

(2) The Levy measure satisfies 

dv{x) = e~ Xi dv[KiX) almost everywhere (4-10) 
meaning that v{B) = f KB e Xi dv(x) for all Borel B. 

(3) The ith coordinate of 7 satisfies 

7i= / x i (l-e^ Xi )dv(x)-\a ii . (4.11) 

Proof. Since 77 is positive integrable, < Ee^ < 00, so that the Esscher 
transform £ l of Q with parameter | and the inverse transform are well 
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defined. According to [32j or Ex. 7.3], £ under £ l has also an infinitely 
divisible distribution, so that 

iff (u) =^v\t( 1 £ \u)-\(u,Au)+ [ (e< u ^-l-t(u,x)l M < 1 )du £ \x) 

for a new vector 7 and Levy measure v . Note that the matrix A is 
invariant under the Esscher transform, see [32j or [34, Ex. 7.3]. 
By Theorem 14. ll fviii) £ G ESDj if and only if 

tpf (u) = (pf {Kju) for all ael". (4. 12) 

By the Levy-Khintchine formula, 

v((Kju) = exp L^, Kju) - ~ (Kju, AKju) 

+ / (e^* T ^> - 1 - i{Kju,x) l^d^ix)) . 

Noticing that (Kju,x) = (u,Kix), changing the variable x to KiX in the 
last integral, using the fQ-invariance of ||| • ||| and the self-inverse property of 
Ki, we see that 

iff (Kju) = exp L{Ka £ \u) - l - {u.K.AKju) 

+ / (e^ -l-%(u,x) I H <i)^ £i (K iX ) ) . 

JR n J 

The uniqueness of the parameters A, u , and 7 of the Levy-Khintchine 
representation for (p £l (see [33l Th. 8.1]) implies that f)4.12p holds if and only 
if 



A — K { AKj , (4.13) 
l £1 =K a £ \ (4.14) 



and the Levy measure v £% is i^-invariant. 

Using the self- inverse property of Ki, representing Ki as the difference of 
the unit matrix and the matrix K[ which has all zeroes apart from the ith 
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column 1, . . . , 1, 2, 1, . . . , 1 with 2 at the ith position, and by equating the 
entries in K[A = A(K-) T , we easily obtain that (14.131) holds if and only if 
a ij = a ji = \ a a f° r ever Y J = j / i. Next, (14.141) holds if and only 

if the ith component of 7 £ * is 0. 

Since the norm ||| • ||| does not change integrability properties in the Levy- 
Khintchine representation, it is possible to replicate the proof from [32] to 
show that the Esscher transform with parameter —\&% leaves A invariant 
while other parts of the Levy triplet are transformed as 

dv{x) = e~* Xi dv e \x) , 

7 = 7 £l + / x{e~* Xi - l)du £ \x) + A( - - e<) . 
■'INI<i 

The latter condition is equivalent to (14. lip , noticing that the ith component 
of A(| &i) is au/2, dv v (x) = e^ Xl dv{x) and ■y £t has zero as its ith component, 
while other components are arbitrary. 
Furthermore, for almost all x, 

dv{x) = e~^ Xi dv s \x) = e~ Xl e~^~ Xi) dv £ \Kix) = e~ Xl dv(Kix) , 

where again we used the fact that the ith component of Kix is — X{ and that 
v £t is Ki- invariant. 

Conversely, the integrability of i] ensures the existence of the Esscher 
transform of Q with parameter | e^. By doing this transform and the converse 
calculations it is easy to verify that Theorem I4.1( viii) applies, i.e. r] = G 
SD 4 . □ 

Since an infinitely divisible random variable £ is symmetric if and only if 
7 vanishes and the Levy measure is symmetric, the above proof is very short 
in the univariate case and immediately yields the corresponding univariate 
result stated in [13 dB], see also [UJ and Corollary 15. 91 

The SD r property of rj implies that the ith component of rj has expectation 
one. If this holds for other components, e.g. if rj forms a martingale, this 
imposes further restrictions on the coordinates of 7, namely 

7j + „ a jj + I ( eXj ~ l ~ x i \A\<i)dv(x) = , j = 1, . . . , n . 

Remark 4.16 (The role of the norm). If we use the Euclidean norm to define 
the truncation in (I4.8p . then this change only affects the value of 7, while 
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A and v remain the same. If 711. 11 denotes the "drift" calculated for the 
Euclidean norm, then 

7||-||=7+ / ^(%||<i-I|*|<i)^> 
so that (14.111) transforms into 

7||-||,i = / x i(%\\x\\<i ~ ^\\x\\<i^ Xi )du(x) - -an . 

Remark 4.17 (Jointly self-dual case). Assume that the conditions of The- 
orem 14.151 hold for each i = 1, . . . , n. The first condition implies that A 
equals up to a constant factor the matrix which has all 1 on diagonal and | 
outside. By applying (14.101) consecutively to coordinates i 7^ j and noticing 
that KiKjKi defines the transposition of the zth and jth coordinates of n- 
dimensional vectors, we see that in this case the Levy measure v is invariant 
under permutations and all components of 7 coincide. 

Remark 4.18 (Finite mean case). Now we also assume that £ has finite mean, 
which is the case if and only if .. >:l ||x||(iz/(a;) < 00, see [531 Cor. 25.8]. Then 
we can rewrite (14. 8B in the following form 

(pf(u)=exp\t(fM,u)-^(u,Au)+ [ (e l{u ' x) - 1 - t(u, x))dv(x)\ (4.15) 

for u G R n , where fi is the Q-expectation of £. Replicating the proof of 
Theorem 14.151 (or by adjusting 7 and using dv(x) = e~ Xi du(Kix)) we obtain 
that £ G ESDj if and only if conditions (1) and (2) of Theorem 14.151 hold, 
while ( 14.111) is replaced by 

IH= I Xi(l -e^ Xi )dv(x) - I- an. (4.16) 
Jr« 2 

Example 4.19 (Log-normal distribution, Black-Scholes setting). Assume that 
rj is log-normal with underlying normal vector £ = log 77, so that 

ipf(u) = exp u) — - (u, Au) j , u G lR n . 

Then 77 G SD.; if and only if the covariance matrix A = (ai m )f m=1 satisfies 
an = au = § a u for I = 1, . . . , n, I ^ i, and fa = -\ an, see (14.161). 
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Finally, rj is jointly self-dual if and only if an = a 2 for all I = 1, . . . , n, 
a im = \° 2 f° r a U ^ 7^ m 5 i- e - f° r cr > the correlations between £j and other 
components of £ are | , and the mean is — ^ for all i = 1, ... , n. The mean 
and covariance matrix of £ are then 



a 2 
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(4.17) 



Remark 4.20 (Square integrable case and covariance). As a consequence of 
Example 14.191 for log-normal rj = G SDj the correlations between the ith 
and other components of £ are ~ y/aa/au (assuming ajj, a^; > 0), while other 
correlations are not affected. In order to relax this correlation structure 
between the ith and other components, it is useful to introduce a jump 
component. For doing that, we assume Jimi^ ||a;|| 2 c/z/(x) < oo, i.e. £ is square- 
integrable. Then the elements of the covariance matrix of £ are given by 



vij = aij + / xiXjdv(x) , 

see [331 Ex. 25.12], i.e. despite of the constrains on the Levy measure given 
in f)4.10p there are various possibilities for the covariance and correlation 
structures. Simple examples can be constructed as in the following remark, 
see also Remark 14.271 and Example 14.311 (for a = 1). 

Remark 4.21 (Levy measures). Assume that £ G ESD« is infinitely divisible 
with the Levy measure v. If v is finite, then the second condition of Theo- 
rem HJ2] means that random vector ( distributed according to the normalised 
v is ESDj itself. In particular, if v is absolutely continuous, its density sat- 
isfies (|4.4|) . An immediate example of v is Gaussian law with the mean and 
variance from Example 14.191 so that is log-normally distributed as in Ex- 
ample 14.191 Since this v is finite, the non-Gaussian part of £ corresponds to 
the compound Poisson law with Gaussian jumps. 



4.4 Quasi-self-dual vectors 

As we have seen, the symmetry properties of random price changes (interpre- 
tation of f] in a risk-neutral case) are considered separately from the forward 
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prices of the assets. In some cases, notably for semi-static hedging of barrier 
options with carrying costs, see [HI Eld]) the symmetry is imposed on price 
changes adjusted with carrying costs a = e A G M n , where usually Xj = r — qj 
for the risk- free interest rate r and qj is the dividend yield of the j'th asset, 

j — 1, ... ,71. 

In view of applications to derivative pricing it is natural to assume that 
all components of rj have expectation one, i.e. Q is a one-period martingale 
measure. Then the random vector 

e x or] = (e Al r7!,...,e A "r7 n ) 

cannot be self-dual with respect to the zth numeraire (resp. for all numeraires) 
unless A, = (resp. all components of A vanish), since the multiplication by 
e A moves the expectation away from one. One can however relate e A o 77 to a 
self-dual random vector by means of a power transformation. 

Definition 4.22. A random vector rj G E n is said to be quasi- self- dual (of 
order a) if there exist AeM" and a^O such that (e A o rj) a is integrable and 
self-dual with respect to the ith numeraire. We then write r\ G QSDj(A, a). 

If rj G QSDj(A,a;), then E(e A, 77i) a = 1 by Lemma [4.81 so that the values 
of a and Aj are closely related to each other. Later in this section, we discuss 
this relation for a special case of quasi-self-dual Levy models. If useful, A 
can also have other interpretations than being the pure carrying costs and 
one can also drop the assumption that 77 is a one-period martingale itself. If 
imposed, the martingale assumption will be explicitly mentioned. 

By Theorem I4.1( iii). 77 G QSDj(A,a) yields that 

E/(e A o V ) = E/(((e A o V ) a )«) = E[/(x,(e A o V )){a lVl ) a ] . (4.18) 

Define random vector £ = A + £, where rj = for £ = (£1, . . . , £ n ). Then 
e A o rj = e^. If we consider the payoff function as a function of asset prices 
St = (Sti, ■ ■ ■ , Stu) with Sxj = So^e^ for j = 1, . . . , n, then (I4.18P can be 
written as 



£>0i , Q q CO Q \\( \ 

■^—{OTl, ■ • • , J T (i-l), OQi, &T(i+l)i ■ ■ ■ 1 ^Tn) j I I 



Fix an asset number % G {1, . . . , n} and assume now that Q is a probability 
measure such that 77 G QSDj(A, a). Since r] a is positive integrable, < 
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Eea^ < oo. Hence, we can define probability measure £ l by 

d£ l dQ 

dQ = Eel ^ ' dl* = E^e-f ^ ' 

i.e. the Esscher transform of Q with parameter ^e^ and the corresponding 
inverse transform. 

It is obvious that 77 G QSDj(A, a) is equivalent to any of the condition of 
Theorem 14.11 for (e A o rj) a = e a< > . The following theorem yields a more direct 
characterisation. 

Theorem 4.23. Let rj a be integrable for some a^O. Then rj G QSDj(A, a) 
is equivalent to one of the following conditions for ( defined from = e x orj = 
e A +«. 

(i) For any payoff function f : E n 1— ► R such that E|/(e^)| < 00 

E/(e c ) = E[f{e Ki <)e aCi ] . (4.19) 

(ii) The distributions of ( and under E % coincide, 
(in) For every u G W 1 , 

<pf(u) = <pf(Kju) 

or, equivalently, 

V ?(u - %e) = ^{kJu - %e)e-^^^ . 

Moreover, if additionally rj is integrable, we have that rj G QSDj(A, a) if and 
only if A4.19\) holds for f being payoffs from basket options with arbitrary 
strikes and weights of assets. 

Note also that all conditions of Theorem 14. 231 can be written conditionally 
on a fixed event or conditionally on a a-algebra, cf. Remark 14.51 The joint 
quasi-self-duality can be achieved by raising the components of rj to different 
powers. 



Proof of Th. \4-23[ For (i) it suffices to note that rj is quasi-self-dual if and 
only if aC, G ESDj and refer to (I4.18P and Theorem I4.1( iii). 

Replace 77 by e A o 77, E[f (^(r]))^) by E[/(x;(e A o rj))(e x *rji) a ], £ i by S\ £ 
by C) and | by | in the proof of the equivalence (iii)-^(vii) in Theorem 14.11 
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to see that (i) is equivalent to (ii). A similar argument yields the equivalence 
of (ii) and <pf '(u) = (pf{Kju) for all u G W 1 as well as the equivalence of 
this equation with 

<p? (u - |* ei ) = cpf (kJu - %e?) (4.20) 

for all u G W 1 . Writing the characteristic functions as Q-expectations and 
using that ( = A + £ yields that 

E exp {x(u - %e h £)+t(u- A) } 

= Eexp {%{Kju - %ei, f ) + t(Kju - A) } . 

Dividing by exp{t{u — f «ej, A)} yields the equivalence of the second statement 
in (hi) and (O0I) . 

If for integrable rj = e^~ x 

E(u + (u,e<)) + = E[(u +(u,x t (e^)) + e a ^} (4.21) 

holds for every (uq,u) G K n+1 we first have that Ee Q ^ = 1 by letting uq = 1 
and Ui — U2 — • ■ ■ — u n — 0. Hence, we can define the measure P by 

so that 

E( Mo + (it, e c » + = E[(«o + («, x i (e^))) + e^] = E P (u + (it, x,(e c ))) + 

for every (w , w) G lR n+1 , i.e., by [2S1 Th. 2.21], under Q and Xj(e^) under 
P share the same distribution. Hence, a payoff function is Q-integrable if 
and only if Ep|/(xj(e^))| < oo and for every Q-integrable payoff-function we 
have 

E/(e<) = E P /(x i (eO) = E[/(x^))e^] , 
i.e. we arrive at (14.191) . The other implication is obvious. □ 

We now use Theorem 14.231 to characterise all quasi-self-dual rj such that 
£ = log// is infinitely divisible with the Levy-Khintchine representation (|4.8|) . 
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Theorem 4.24. Let the random vector £ = log 77 be infinitely divisible under 
Q with the generating triplet (A, is, 7) and let n a be integrable for some a ^ 0. 
Then 77 G QSDj(A, a) if an only if the following three conditions hold. 

(1 ) The matrix A = (%-)y=i satisfies = = | ctjj for all j = 1, . . . , n, 
3 7 1 i- 

(2) The Levy measure satisfies 

dv(x) = e~ aXi dis(Kix) almost everywhere (4.22) 
meaning that is(B) = f KB e aXi dis(x) for all Borel B. 

(3) The ith coordinate of 7 satisfies 

7i= / Xi(l - e^ Xl ) dis(x) - - Xi . (4.23) 

Proof. Denote ( = X + £. Since < Eef ^ < 00, the Esscher transform E l 
of Q with parameter ~ej and the inverse transform are well defined. There- 
fore, ( under 8 l has also an infinitely divisible distribution. By using Theo- 
rem H723]^iii) instead of Theorem I4.1( viii) and replacing £ % by £ % , £ by £, \ 
by f (14.81) in the proof of Theorem I4.15[ we obtain (1), (2), and 

H = I Xi(l - e% Xi )dis(x) - \n 
J\M\<i 1 

for the generating triplet of ( under Q. Since £ = ( - A we only have to 
adjust 7j by — Aj to finish the proof of the first implication. 

The integrability of i] a implies the existence of the Esscher transform of Q 
with parameter ^e^. By doing this transform and the converse calculations it 
is easy to verify that Theorem I4.23( iii) applies, i.e. 77 = G QSDj(A, a). □ 

Note that condition (1) is identical to Theorem 14.151 (1) . As a consequence 
of Theorem 14.15( 2) and I4.24T 2). we immediately get the following results. 

Corollary 4.25. Let the random vector £ = logn be infinitely divisible under 
Q with non-vanishing Levy-measure is. Then 77 cannot be quasi-self-dual of 
two different orders with respect to the same numeraire. 
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Corollary 4.26. If £ t , t > 0, is the Levy process with generating triplet 
(A, u, 7) that satisfies the conditions of Theorem \4.24l then G QSDj(At, a) 
for all t > 0. 

Proof. It suffices to note that ff t {u) = for all t > and raise the 

corresponding identity from Theorem I4.23( iii) into power t. □ 

Remark 4.27 (Levy measures in the quasi-self-dual case). In order to con- 
struct a Levy measure v satisfying (14.221) . note that 

e^ Xi dv[x) = e^^ KiX ^dv(Kix) , 

meaning that the measure u with density ^j-(x) = e% Xi is fQ-invariant. 
Therefore, in the background one always needs to have a .^-invariant Levy 
measure. 

Since the Lebesgue measure on K n is fQ-invariant, a simple example of 
z/o is provided by the Lebesgue measure restricted onto Br, where Br = {x : 
IMI — R} is the ball of radius R in the |||-||| norm. A further implication of the 
i^j-invariance property of the Lebesgue measure on R n is that the Lebesgue 
density p uo of an absolutely continuous i^-invariant measure v$ is also K r 
invariant, i.e. p vo {x) = p U0 (Kix) for almost every x G M n . Then (14.221) can 
be equivalently written as p u (x) = e~~ aXi p v (Kix). Clearly, condition (14. 9 p is 
always satisfied for a finite v without atom at the origin, which then yields 
the compound Poisson part of £ from 77 = G QSDj(A, a). The integrability 
condition on n a additionally requires that 

/ e ax i e -f x * dv Q {x) < 00 , j = l,...,n, 

J\\x\\>l 

see [33, Th. 25.17]. 

Remark 4.28 (Determining a from the carrying costs in the risk-neutral case) . 
Assume that Et^- = 1 for all j = 1, . . . , n and r\ G QSDj(A, a) with given A. 
Since ip^(-tej) = Enj = 1, we see that 

7j ■• = ~ I {e Xj -1-Xj l M \<i)dv(x) - - a j3 , j = l,...,n. (4.24) 
If a = 1, then the above condition for j = i yields (14.111) (or (14.231) for a = 1 
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and A = 0). Indeed, it suffices to check that 
/ (1 - e Xi + x^ Xi %|<i) du(x) 

JR n 

= J (l-e^+Xie^' I| x |<i)e" x< du(K iX )+ J (l-e^+^e^ 1 I| x |<i) du(x) 

{xi<0} {xi>0} 

= J {e yi -l-y i e^ yi l M < 1 )dp{y)+ J {l-e Xi + Xi e^ Xi l lxl <t) dv{x) = . 

{Vi>0} {xi>0} 

However, for non-vanishing A we need to combine (I4.24p with (14.231) to see 
that a must satisfy 

f \ f a Ct 

-I (e Xi -l-x i li x i< 1 )du(x)--au= / - e* Xi ) dv{x) - -a u - A, , 
jR n 2 J I 



INI<i 

or, equivalently, 



a u a = a {i -2X { + 2 / (e a ' 1 - 1 - x^ 1 1\ x \<x)<k/(x) 

JR n 



(4.25) 



It should be noted that in the Levy processes setting from Corollary 14.261 the 
values of a calculated for all t > coincide. 

Remark 4.29 (Finite mean case). Assume that 77 = G QSDj(A,«). If, as 
in Remark 14.181 £ has finite mean, then (14.231) is replaced by 

Mi= / Xi(l - e% Xi ) dv(x) - % u - Xi, (4.26) 

JR n ^ 

where /1 is the expectation of £. If Q is a martingale measure for rji, then 
*ej) = Ee fi = 1 yields that 

IH = - I (e Xl - 1 - Xi)dv{x) -\a u . (4.27) 
JR n 2 

Combining (Q5]) with (Q7|l yields 

a^a = a« - 2Ai + 2 / [e Xl - 1 - a^et^cfe^) 

= an - 2\i + 2 / (e^ - 1 - Xie^du^Xi) , (4.28) 
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where Vi is the marginal Levy measure defined by Vi{B) = u({x el": Xj G 
B}) for Borel B C I, ^ 5, see [331 Prop. 11.10]. 

Compared to (I4.25p . Equation (14.281) yields a considerable simplification 
in calculating a. Since z/, is the Levy measure corresponding to r/i, it is 
possible to calculate a from only the distribution of the zth component of rj 
and the corresponding carrying costs Aj. 

In the purely non-Gaussian case (i.e. if A vanishes) it is useful to write 
the integral in (14.281) as its principal value. Then the principal value of the 
integral of XiC^ Xi vanishes, since dvi{xi) = e~^ Xi dvQi{xi) for a symmetric 
measure vqh, and 



If i>Qi has a finite Laplace transform tp on the real line, then a solves 



Example 4.30 (Log-normal model with carrying costs). By Corollary I4.25[ 
among all log-infinitely divisible distributions only the log-normal one can 
be quasi-self-dual of two orders with respect to the same numeraire. Apply- 
ing (I4.25P for the univariate log-normal case with an = a 2 > (and vanishing 
v) yields that 

a = l , 

a 1 

as stated in [8j [9], [11]. Hence, the univariate log-normal distribution in the 
Black-Scholes setting is self-dual and quasi-self-dual of order 1 — % at the 
same time. By (I4.25p . this is also true for multivariate log-normal models 
from Example 14. 191 being self-dual with respect to the ith numeraire, i.e. this 
distribution is at the same time quasi-self-dual of order a = 1 — 2Xi/au with 
respect to the same numeraire. 

Example 4.31 (Determining a for non-trivial Levy measures). Start with the 
univariate case (i.e. n = 1) and choose v from Remark 14.271 to be the centred 
Gaussian measure with variance (3 2 > 0. If normalised to have the total mass 
one, v becomes the density of the normal law with mean — and variance 



Ai = f(e Xi - l)dvi{xi) = 4{e Xi - l)e * Xi dv 0i {xi) 
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(3 2 . Solving (14.251) or equivalently (I4.28P for this particular measure v and 
an = a 2 > yields that 

a = JL ( 2 Lambertw(^exp {^±11})^ + /?V - 2/3 2 A - 2/3 2 ) , 

where Lambert W(x) = g(x) is the principal branch of the Lambert W func- 
tion that satisfies g(x)e 9 ^ = x for all x. In the purely non-Gaussian case 
the required power is given by 

« = l-^log(l + A). 

In the multivariate case we start with Uq being the centred Gaussian law 
having positive definite covariance matrix B that satisfies Theorem 14.24( 1) 
for some fixed i and define measure v with density 

dv _* x . 
-— (x) = e 2 8 . 

Then (14.221) holds and the zth marginal z/j of v has the density e~^ Xi with 
respect to the ith marginal of uq, the latter being the centred Gaussian law 
with variance f3 2 = b„. Since the ith marginal for the normalised v coincides 
with the Levy measure constructed above in the univariate case, we obtain 
the same a as in the univariate case with (3 = \fbii. 



5 Distributions of self-dual random variables 
5.1 Characterisation and examples 

In this section we specialise the results from Section H] for studying self-dual 
random variables. Denote by F(x) = P(i] > x) the tail of the cumulative 
distribution function of a positive random variable r\ and by 

F^z) = [ F(t)dt, z>0, 
Jo 

the integrated tail. Note that -Fj(O) = and -Fj(oo) = 1 in case Er] = 1. 
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Theorem 5.1. An integrable positive random variable 77 is self-dual if and 
only if Fj (00) = 1 and 

zF^z' 1 ) = F^z) for allz>0. (5.1) 

Proof. It is easy to check that 

Fi(z) = Emin(?7, z) , z > . 

Now apply Theorem 12.41 or Theorem 14.11 (iii) to the payoff function f(rf) = 
min(r), z) to see that 

Fj(z) = Emin(?7, z) = E[min(r/" 1 , z)rj\ = Emin(l, 277) = zFi(z~ l ) . 

In the opposite direction, (15. II) yields that 

E max(r/, z) = E[r/ + z — min(?7, z)] = E[l + zrj — min(l, zrj)] = E max(l, zrj) , 

i.e. by rescaling (cf. Remark I4T21) we arrive at the self-duality property (13. 4p . 

□ 

Theorem 14.71 in the univariate case yields the following result, which is 
known from [371 Ex. 8]. 

Corollary 5.2. Let n be a positive integrable random variable with distri- 
bution Q. 

(a) Iff] is absolutely continuous with probability density p v , then rj is self- 
dual if and only if 

p v (x) = x _3 ^(x _1 ) for almost all x > . (5.2) 

If £ — log 77, the self-duality of 77 (i.e. the exponential self-duality of £) 
is equivalent to 

p%(x) = e~ x p^(—x) for almost all x G M . 

(b) Iff] has a discrete distribution, then 77 is self-dual if and only if Q(ry = 
a; -1 ) = xQ(r] = x) for each atom x of rj. 
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Clearly, if the density p v is continuous, then (15.21) holds for all x > 0. For 
instance, the probability density of the log-normal distribution of mean one 
satisfies ( 15. 2j) . It is also satisfied by mixtures of log- normal densities that 
appear in the (uncorrelated) Hull- White stochastic volatility model, see [211 
Th. 3.1]. The self-duality property of stochastic volatility models is explored 



Example 5.3 (Log-normal model). If St = Ft] has the log-normal distribu- 
tion, the Black-Scholes formula yields that 



able. Note that the conventional Black-Scholes formula is obtained by sub- 
tracting k from (15.31) and then discounting. By looking at the right-hand 
side of (15.31) it is easy to see that it is symmetric with respect to F and k, 
i.e. rj is a self-dual random variable. 

The right-hand side of (15.31) defines a (symmetric) norm on R 2 ^ called 
the Hiisler-Reiss norm of x = (k,F), see [27]. Thus, the derivative given 
by the maximum of the asset price and the strike has the price given by 
the discounted norm of the vector composed of the forward and the strike. 
Notably, expression (15.31) appears in the literature on extreme values, see [22J , 
as the limit distribution of coordinatewise maxima for triangular arrays of 
bivariate Gaussian vectors with correlation g(n) that approaches one with 
rate (1 — g(n)) logn — > A 2 G [0, oo] as n — > oo. 

In order to construct further examples of probability density functions 
p r] that satisfy (I5.2p it suffices to define p v (x) for x > 1 and then extend it 
for x G (0, 1) using (15.21) with a subsequent normalisation to ensure that the 
total mass is one. Clearly, one has to bear in mind that E77 = 1 presumes 
the integrability of xp v (x) (alongside with p n {x) itself) at zero and infinity. 

Example 5.4 (Self-dual random variables with heavy tails). The log-normal 
distribution has a light tail at infinity. It is possible to construct a self-dual 
heavy-tail distribution by setting 



in [TH Th. 3.1]. 



Emax(F??, k) = F<$>(d') + k$(d") , (5.3) 

where k, F > 0, 





(5.4) 
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for 7 > —1, where c 7 = (1 + 7) (2 + 7)/(3 + 27) normalises the probability 
density. 

Example 5.5 (Discrete self-dual random variable). If i] takes values |,1,2 
with probabilities ~, | , I, then Corollary 15.2( b) implies that n is self-dual. 

Remark 5.6. If 77 is not self-dual, then (13.41) is clearly violated, but the result- 
ing inequalities can not be the same way around for every k, F > 0. Without 
loss of generality assume that F — 1. Then E max(^, )}) < Emax(^, 1) for 
all k > with the strict inequality for some k = k leads to a contradiction, 
since 

Emax(/co,?7) < Emax(/co?7, 1) = A;oEmax(?7, A;^ 1 ) 

< A; E max(A; _1 r], 1) = Emax(j), /cq) . 



5.2 Moments of self-dual random variables 

It is immediate that all self-dual random variables have expectation one. 
Carr and Lee [Til Cor. 2.6] show that 

Erf = ErT n+1 , n > 1 . (5.5) 

In particular, if En 2 < 00, then 

Cov(t7~\ 77) = 1 - E/T 1 = ( E? ?) 2 - E ^ 2 = -Var(r/) . 

Theorem 5.7. Each non-trivial self-dual variable 77 with finite third moment 
has a positive skewness E(?7 — En) 3 / (Var(n)) 3 ^ 2 . 

Proof. In view of (15.51) . 

E(r? - En) 3 = E(n 3 - 3En 2 + 2) 

= E(n 3 - 3Er] 2 + 2 + 6(77 - 1) + r]- 1 - n 2 ) 

{n-lfirt + r,- 1 - 



E 



> 0. 



This also shows that the skewness vanishes if and only if r\ — 1 almost 
surelyQ □ 



1 The authors thank a referee for suggesting the current proof. 
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Remark 5.8 (Product of self-dual variables). If rjx and 77 2 are two independent 
self-dual random variables, then 

E max (/c, 771 7/2) = E[E(max( k, 771772) |»7i)] 

= E[E(max(/c77 2 , ^1)^2)] = Emax(fc77i772, 1) , 

i.e. the product 771772 is self-dual. By taking successive products it is possible 
to construct a sequence of self-dual random variables, whose logarithms build 
a random walk. Note however that the values of this random walk at different 
time points are not jointly self-dual, cf. Remark 14.121 



5.3 Exponentially self-dual variables 

Theorem 14. lT viii) implies that £ is exponentially self-dual if and only if the 
characteristic function (p® satisfies 

<pf{u--t) = (pf(-u--t), ueR. 

If £ has an absolutely continuous distribution, Corollary 15.2( a) yields that £ 
is self-dual if and only if e^ y p^{y) is an even function of y. 

If £ is also infinitely divisible, then its distribution is characterised by 
the Levy triplet (<r 2 ,z/, 7). Note that in the univariate case K\x = —x, the 
norm (14. 7p becomes the Euclidean one and A reduces to a single number a 2 . 
Theorem 14 . 1 51 yields the following univariate result, known from Fajardo and 
Mordecki [HJ[T8]; to see that their "drift" with truncation function l|a;i<i is 
equal to 7 from Corollary 15.91 use e~ x dv{— x) = du(x). The latter condition 
on the Levy measure appears also in Carr and Lee [TTJ Th. 4.1]. 

Corollary 5.9. An integrable random variable rj = with £ being infinitely 
divisible represented by the Levy triplet (& 2 ,v, 7) is exponentially self-dual 
if and only if dv(x) = e~ x du(—x) and 

r 2 

7 = / x(l - ea x )du(x) . (5.6) 

J\x\<i 2 

If £ is integrable, then ( I5.6'|) can be replaced by the following condition on its 
expectation 

fJi — I x(l - e^ x )du(x) — — . 
Jr 2 
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While Corollary 15.91 is obtained as a univariate version of Theorem 14.151 
it is alternatively possible first to describe the (univariate) dual market in 
terms of its generating triplet, and then ensure that the generating triplets of 
the original and dual markets coincide, implying that rj is self-dual, see [17]. 
The latter approach also describes the dynamics of the dual market in the 
univariate case. 

If E|£| 3 < oo, then [EH Prop. 3.13] yields that 

p POO 

E(£-E£) 3 = / x 3 dv(x) = / x 3 (l-e x )du(x). 
Jr Jo 

Thus, the skewness of exponentially self-dual £ is negative except in the log- 
normal case, where it is zero. 

5.4 Quasi-self-dual variables and asymmetry correc- 
tions 

Let St = S ar] for S , a > with r\ being a general positive random variable, 
so that the forward price is given by F = S^a. Assume that rj is absolutely 
continuous with non- vanishing density p v and Erj = 1. Then it is possible to 
find a function q ar) such that 

Ef(S T ) = E[f(So/(a V ))q av (ari)} = E[/(F/(a 2 r / ))g a ,(a7 ? )] 

= E[f((So) 2 /ST)qa V (ar ] )] (5.7) 

for each function / : IR + i— > R such that /(St) is integrable. Indeed, it 
suffices to choose 

, s Pa-nix' 1 ) Ps T ( x ' ls o) 

Variix) " 



X 2 p av (x) X 2 ps T (xS ) ' 

By choosing x = arj — St/S we arrive at 

(S )\,St,-2 PsA(So) 2 /S t ) 



E/OSr) = E 



St Sq Ps t \St) 

Apart from trivial cases, the density ps T of St depends on T. In view of 
applications to semi-static hedging described in [TT] it is beneficial if the 
correcting expression 

p St ((S ) 2 /x) 

QsAx) = ^ 

Ps t {x) 
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at any time t G [0,T] depends only on x and S but not on t. This is the 
case if rj is self-dual with no carrying costs (then qs t ( x ) — ( x /So) 3 , x > 0, by 
Theorem \2A\i , or quasi- self- dual with parameters a = e x and some a / 0, 
being the case if and only if qs t (x) = (x/ So) 2+a , x > 0. In the latter case (15.71) 
turns into 

Ef(F V ) = E\f{^ r )a a V c 

By letting f(x) = (x—k)® and noticing that E(ar]) a = 1 in the quasi-self-dual 
case, this implies the following property 

E(Fi] - k)% = a~ a E{F - ka 2 r])% = E V a E(F - k V (Er] a )~ 2/a ) + , 

which can be termed as the power put-call symmetry. 



6 Barrier options and semi-static hedging 

6.1 Time-dependent framework 

Consider a finite horizon model with the asset prices given by 

S t = S o e tx o Vt = S o e tx+ ^ = (S 01 e tx ^\ S Qn e tXnHt ") , t G [0, T] , 

where A G M n represent deterministic carrying costs and all components 
of r]t = e& are martingales with t G [0, T], being a Levy process. Fix 
i G {l,...,n} and assume that r\ t G QSDj(tA,a) for every t G [0,T]. This 
condition is satisfied (with a = 1 and A = 0) for all exponentially self-dual 
Levy models with no carrying costs analysed in Section H~3l and for quasi-self- 
dual Levy models from Section fl~4l for non- vanishing A, see Corollary 14.261 

Let t be a stopping time with values in [0, T] and let 3v be the correspond- 
ing stopping cr-algebra. Since £ 4 is a Levy process, (£ t ,£t) an d ((,t,£,t + £,'t- t ) 
share the same distribution, where ^, t G [0,T], is an independent copy of 
the process * e [0,T]. Hence, (S T ,S T ) and (S T ,S T o e A(T " T)+ ^- T ) also 
coincide in distribution. Then 

E[/0Sr)|&] = E[/(5r)|5 r ] = E[/(5 T o e A ^ + ^)|^] 

= E[/(S T oe A < r ^+&-)|3v], 
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where / is any integrable payoff function. The quasi-self-duality of r]' T _ T = 
e^ T - T with respect to the ith numeraire adjusted for conditional expectations 
(see Remark [4. 5 p yields that 

E[f(S T )\$ T ] = E[f(S T o e ^WM + ^) )e ^(M + ^ r)i ) |5T ] ; 

whence 



= E 



. STlS T i ST(i-l)S T i ST(i+l)S T i STnS T i\ / Sri , ^ 

J \ ri ' ' ' * ' c ' ~0 ' C ' ' ' ' ' C / I ~C ' 



(6.1) 



cf. Remark 14 .41 If S T i = H almost surely for a constant H, then 

nf(s T )\dr] 



E 



St\H St(i-1)H H 2 ST(i+l)H St,iH\ / STi 



J \ q > • • • ) q ) "7< i q ) • • • ) q I I rr / 

>->Ti >->Ti >->Ti iJTi >->Ti / ^ n 



(6.2) 



Identity (16. 2p for a = 1, A = yields the self-dual case, and in the univariate 
case n = 1 corresponds to [TTj Eq. (5.3)]. Classical examples with trivial 
carrying costs (i.e. A = 0) are options on futures or options on shares with 
dividend-yield being equal to the risk-free interest rate. For the univariate 
quasi- self- dual case, see [HI Cor. 5.10]. 

Remark 6.1. Instead of the self-duality property, it is possible to impose (16.11) 
for stopping times r G [0,T] that might appear in relation to hedging of 
particular barrier options. This observation leads to extensions for indepen- 
dently time-changed multidimensional Levy processes by means of condition- 
ing arguments described in [Til Th. 4.2, 5.4]. Further models without jumps 
can be obtained on the basis of the multivariate Black-Scholes model with 
characteristics described in Example 14.191 by applying independent common 
stochastic clocks being continuous with respect to calendar time. 



6.2 Multivariate hedging with a single univariate bar- 
rier 

Assume a risk-neutral setting for a price process S t , t G [0,T], and fix a 
barrier level at H > such that So% 7^ H with given i G {1, . . . ,n}. For 
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simplicity of notation, define function Xj : E n t— > E ra acting as 

ki{ST, H) = — — (Sti, • • • , Sr(i-i), H, Stu+i), • • • , Sm) ■ 

Define E ti to be the (closed) line segment with end-points and Su and let 

t h = mf{t >0:He E ti } and \ = ^t h <t , 

cf. [TT| Sec. 5.2] who used a bit different way to handle the two cases when the 
initial price is respectively lower and higher than the barrier. Furthermore, 
assume that the asset price dynamics satisfy (16.11) for the stopping time r = 
th and that S TB i = H a.s. on the event that {r# < T}, what is guaranteed, 
e.g. by the sample path continuity of the ith component of S t , t > 0. In 
case of discontinuous Levy processes, the symmetry condition f)4.22p on the 
Levy measure implies the presence of jumps of both signs, so it is much more 
difficult to ensure that S TH = H a.s. 

Take any integrable payoff function / and consider an option with payoff 
xfi^r), the knock-in option with barrier H for the zth asset. In order 
to replicate this option using only options that depend on the terminal value 
St consider a European claim on 

f(S T ) t H&n+ (^Y f^S T ,H)){t HeBTi - l STi=H ) . (6.3) 

Here one has to bear in mind that this is only practicable provided that 
the considered claims are liquid or can be replicated by liquid instruments. 
However, there is a fast growing literature about sub- and super-replication 
of multiasset instruments, see e.g. [23] and the literature cited therein. 

On the event that {th > T}, the claim in (I6.3P expires worthless as 
desired. If the barrier knocks in, we can exchange (16.31) for a claim on /(St) 
at zero costs. To confirm this, define E# to be the (closed) line segment with 
end-points S 0i and H 2 S^ . Note that H ^ E ti if and only if H e E ti \ {Su}. 
Hence, on the event that {r# < T}, by (16.21) . we have 

E[/(St)|3v] = W(St) 1W Ti |3v] + W(St) I^b„ ISA 

= B[f(S T ) l H ez Tl \dr] + E\(^) a f(Sc i (S T , H)) l H ^ n \!g T 



+ E 



(^) a f(Ki(S T ,H))(l HeETi - Ks Ti =H)\$r 
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For simplicity we assume from now on that Sti has a non-atomic distri- 
bution, so that (16.31) becomes 



(f(Sr) + (-§ ) a f(HS T ,H))) l H&Ti . (6.4) 

Consider general basket call /(St) = (J2 7 j=i u jSTj — k)+. By (16. 4ft the 
hedge for the knock-in basket call with payoff function xf(Sr) is given by 
the derivative with payoff function 

3=1 3=1, 3'¥* 

which depends only on St- 

If (16.11) holds with a — 1, this hedge becomes 

n , n 

{(j2 u J S Tj ~ k ) + + ( U * H ~ (ff STi ~ Yl u i ST i))J ^H& Ti , (6.5) 

3=1 3=1, j¥* 

being the sum of a basket call and a spread put with knocking condition 
depending only on the ith component St% at maturity. 

In some cases the knocking condition at maturity can be incorporated 
into the payoff function. For this, note that we can write all integrable 
payoff functions in the form 

f(S T ) = fo(S T ) ls T ee , where 9 = {x : f(x) + 0} , 

with 6 possibly being E n . For example, the basket call (J2 1 j=i u j^Tj — k) + 
can be written as (J2 7 j=i u jSTj — k) ^Y^=i u jS T j>k- Of course if H En 
would imply that ki(ST,H) ^ 9 and St ^ O at the same time, then it is 
possible to omit I# g ~ Ti in (16. 4p . but this is not the case for standard payoff 
functions. If H e En implies that St ^ Q (resp. ki(Sx,H) 0) then 
the first (second) summand in (16. 4p is always zero. If furthermore H ^ En 
implies that ki(ST,H) 6 (St ^ ©) then we can omit the first (second) 
summand in (16. 4p and hedge with the second (first) summand without the 
knocking condition I# g = T ., i.e. in (16.51) we can hedge with a conventional 
basket option. 

Example 6.2. Consider a bivariate price process (S t i,S t 2) in a risk-neutral 
setting satisfying (16.11) with a — 1 and i — 1 for the stopping time r = t h = 
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inf{£ : S t i < H} with barrier H such that < H < Soi- First assume again 
that Sti can not jump over H . For the spread option 

/(Sti, S T2 ) = (aS T i - bS T 2 - k)+ , a, b > , 

assume additionally that aH < k and define x — ^-t h <t- By using the hedg- 
ing strategy described in (16.41) and aH < k we obtain a henge for xf(Sri, <St2) 
by 



(aS T i - bS T 2 - k) + l He[STuSoi] + (aH - jjS T i ~ bS T 2) ^He 



[S T i,Soi] 



aH — —S T1 — bS T ^j 1_h- 6 [5 t1i 5 01 ] = [aH — —S T1 — bS T2 

i.e. it is possible to hedge with a basket put. Therefore, the related knock- 
out option can be hedged with a long position in the spread call with payoff 
function /(Sti, SV2) — (ciSti — bSx2 — k) + and a short position in the above 
hedge. Note that we only assumed that b > so that the knock-in level can 
but need not be deep out-of-the-money. If S t i can jump over the barrier H we 
get a super-replication in case of the knock-in option and a more problematic 
sub-replication in case of the knock-out option. 

Assuming (16. ip for the stopping time th with q^I, where S t \ does not 
jump over H, the hedge for the knock-in option has to be modified as 

STi\ a ~ 1 ( k 

aH — —Sti — bST2 



H ) \ H 

while the modification for the knock-out is now obvious. This example can 
easily be extended in a higher dimensional setting as long as all risky assets 
without knocking barrier enter the payoff function with a minus sign. 

Example 6.3. Consider a bivariate price process in a risk-neutral setting. 

Assume that S01 > k > and that (16. ip holds for r = = inf{t : Sti < k}, 

a = 1, and i — 1, while Sti can not jump over k. Introduce (possibly 
negative) payoff function 

g(S T i, S T 2) = (S T i - k) + (S T 2 A (Sri - k)) , 

where a A b = min(a, b). 

By (16.41) with a = 1 we obtain a hedge for lL Tk < T 9(Sti, Sy 2 ) as 

\{(S T i-k) + (ST2A(STi-k))) + {(k-STi) + (ST2A(k-S T i)))}l k e 
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Here we can get rid of the indicator function by noticing that the above 
payoff function can be written as 

((S T i - k) + (S T2 A (S T i - k))) - (S T i + S T2 - k) + + (k + S T2 - S T1 ) + . 

Furthermore, for the related knock-out option we get a hedge given by a long 
position in the basket call with payoff function (Sti + St2 — k) + and a short 
position in the put spread with payoff function {k — (Sti — St2))+- 

6.3 Examples of hedging in jointly self-dual cases 

In this section we create hedges for more complex instruments and bivariate 
models satisfying (16. ip for two different numeraires, e.g. for jointly self-dual 
exponential Levy models with generating triplets satisfying the conditions in 
Remark H~TT1 

Example 6.4 (Options with knocking-conditions depending on two assets). 
Assume a risk-neutral setting for a price process S t = (S t i,S t2 ), t G [0, T], 
where (16 .ip holds for both assets with a = 1 and the subsequently defined 
stopping times. Furthermore, let k x , k y > be constants such that k x < 
Soi,So2 < k y and define the stopping times = inf{t > : Su < k x }, 
Tiy = inf{£ > : Su > k y } and the corresponding stopping a-algebras $ Tix , 
$ Tiy , i = 1,2, as well as the stopping time r = t 1x At 2x with the corresponding 
stopping cr-algebra $ T . Assume also that the price processes can not jump 
over the barriers k x and k y respectively. 
Consider the claims 

X = (Sti — St2 — k x ) + ~& t=T1x <t +(St2 — Sti — k x ) + 1 t1x ^t=t 2x <t , 
Y = (k y - S T i - St2)+{% 1v <t<t 2v - ^T 2y <T<r ly ) , 

i.e. A is a knock-in spread option on the difference between the share which 
first hits k x and the other one only being knocked-in if at least one share hits 
k x before T. At maturity, Y is a long position in a basket put if and only 
if the first but not the second asset price hits the price level k y and a short 
position in the same basket put if and only the second but not the first asset 
hits k y before T. 

The claim X can be hedged with a long position in the European basket 
put with payoff (k x — Sti — St2)+- The claim Y can be hedged by entering a 
long position in the spread option with payoff (SVi — SV2 — k y ) + along with 
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a short position in the spread option with payoff (St2 — Sti ~ k y ) + . To see 
that, apply identity (16. 2p for a = 1, so that 

n(k z - S T1 - S T2 ) + \$ Tlz ] = E[(S T1 - S T2 - k z ) + \$ Tlz ] , (6.6) 
n(k z - Sti ~ Si^+I&s,] = E[(5 T2 - S T1 - k z ) + \$^} , (6.7) 

z — x,y, while the value of the European spread option with payoff (Sti — 
St2 — k) + (resp. (St2 — Sti — fc)+) remains unchanged by applying (16.21) at 
t 2 (resp. n). 

As far as X is concerned we have that in case where {r > T} neither X 
is knocked in nor the basket-put in the hedge portfolio is in the money, since 
Sti, St2 > k x . If {r = r lx < T}, by (16.61) we can exchange the long position 
in the hedge portfolio for the needed spread, if r lx ^ t = r 2x < T the same 
is true due to (16.71) . 

As far as Y is concerned we have that if {ji V > T}, i = 1,2, both 
instruments in the hedging portfolio are out of the money since Sti, St2 < k y , 
i.e. have payoff zero as Y. On the event that we first have {ri y < T} we can 
change the long position in the spread option with payoff (Sti — St2 — k y ) + in 
a long position in the basket put with payoff (k y — Sti ~ 5r 2 )+ while letting 
the short position unchanged. Provided that additionally {j2 y G [t^jT]}, 
by (16. 7p . we can also exchange the short position for the same basket put, so 
that we can close our positions as required, otherwise, i.e. {r 2y > T}, the long 
position in the hedge portfolio yields a potentially needed payoff of Y while 
the short position still matures worthless. On the event that we first have 
\j2y < T}, by (16. 7p . we can change the short position in the needed basket 
put while letting the long position unchanged. If furthermore {ri y G [r 2y , T]} 
we can again close our position due to (16.61) . Otherwise, unlike the long 
position, the short position in the hedge portfolio may be in the money at 
maturity but in that case Y at maturity would be a long position for the 
hedger with the same payoff. 

Assuming that (16.11) holds for ot\ resp. a 2 with respect to the correspond- 
ing components (where other assumptions remain unchanged), we can ap- 
ply (I6.2p in the same way to see that the hedge should theoretically be mod- 
ified to a long position in the European derivative with payoff (Sti — St2 ~ 
k y ) + (k~ 1 5Vi) ai_1 along with a short position in the European derivative with 
payoff (S T 2 - Sti - k y ) + (k~ l S^)^ 1 ■ 

For creating semi-static hedges of barrier spread-options with certain 
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knocking conditions, e.g. claims of the form 

Z = (S T1 — S T 2 — k) + Is tl >St 2 ,v*e[o,T] > k > , 

in equal carrying cost cases the full strength of the joint self-duality is not 
needed. It suffices to assume exchangeability being implied by the joint 
self-duality, see Corollary 14 . 1 01 and [28J for details including model character- 
isations, weakening of the exchangeability assumption and hedges for several 
related derivatives. 



6.4 Semi-static super-hedges of basket options 

The following super-hedges may be quite expensive for replication purpose. 
Thus, they seem to be more useful if one would like to speculate with a basket 
option and get some additional money by writing a different knock-in basket 
option, where the maximum loss should be limited to the initially invested 
capital. 

In the sequel we work in the same setting as in Section 16.21 with the 
additional assumptions that a = 1 and Soi > H > 0. Define again the 
stopping time r# = inf{t : Su < H} and let $ TH be the corresponding 
cr-algebra. Consider the knock-in basket option with the following payoff 
function 



X 



n 

( u o^ T i ~ , k, Ui > 0, Uj £ M for j — 1, . . . , n , j ^ i . 



where x = ^t h <t- 

By (16.21) for a = 1 we have 



E[(53«i5ri-A:) + |^]=E[( £ Uj S Tj - ^S Ti + « 4 H) + |ff^ 

Hence, the maximum loss of buying the basket option in the right-hand side 
of the above equation and short-selling the initial knock-in basket call does 
not exceed the initial costs of this strategy. Note that in this setting, jumps 
over the barrier H would add a further aspect of super-hedging. 

If (16.11) holds for a ^ 1 where S ti can not jump over H, use again (16. 2p . 
The long position in the described strategy would then be determined by the 
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payoff function 

(if) ( ZJ u j S Tj -^S Ti + u i H) + , 

j •: J'/' 

while the short position remains unchanged. 
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